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ABSTRACT 

Context. With the current study we aim at understanding the properties of radio emission and the assumed jet from four radio-loud and 
y-ray-loud narrow-line Seyfert 1 galaxies. These are Seyfert 1 galaxies with emission lines at the low end of the FWHM distribution. 
Aims. The ultimate goal is twofold: first we investigate whether a relativistic jet is operating at the source producing the radio output, 
and second, we quantify the jet characteristics to understand possible similarities with and differences from the jets found in typical 
blazars. 

Methods. We relied on the most systematic monitoring of radio-loud and y-ray-detected narrow-line Seyfert 1 galaxies in the cm 
and mm radio bands conducted with the Effelsberg 100 m and IRAM 30 m telescopes. It covers the longest time-baselines and the 
most radio frequencies to date. This dataset of multi-wavelength, long-term radio light-curves was analysed from several perspectives. 
We developed a novel algorithm to extract sensible variability parameters (mainly amplitudes and time scales) that were then used 
to compute variability brightness temperatures and the corresponding Doppler factors. The jet powers were computed from the light 
curves to estimate the energy output and compare it with that of typical blazars. The dynamics of radio spectra energy distributions 
were examined to understand the mechanism causing the variability. 

Results. The length of the available light curves for three of the four sources in the sample allowed a firm understanding of the 
general behaviour of the sources. They all display intensive variability that appears to be occurring at a pace rather faster than what 
is commonly seen in blazars. The flaring events become progressively more prominent as the frequency increases and show intensive 
spectral evolution that is indicative of shock evolution. The variability brightness temperatures and the associated Doppler factors are 
moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case 
clearly implies a quiescent jet underlying recursive flaring activity. Finally, in one case, the sudden disappearance of a y-ray flare 
below some critical frequency in our band needs a more detailed investigation of the possible mechanism causing the evolution of 
broadband events. 

Conclusions. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are 
powered by relativistic jets, such as intensive variability, spectral evolution across the different bands following evolutionary paths 
explained by travelling shocks, and Doppler factors indicating mildly relativistic speeds. 
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1. Introduction 

Seyfert galaxies were first_jdentificd as a distinct class of ex- 
tragalactic systems by [Seyfert! (1 19431) . who studied the nuclear 
emission of six “extragalactic nebulae” with “high excitation nu¬ 
clear emission lines superposed on a normal G-type spectrum”. 
The lines appeared to be broadened, reaching widths of about 
8500 km s^ 1 . [Seyfert! also noted that the maximum width of the 
Balmer emission lines increased with the luminosity of the nu¬ 
cleus as well as the ratio between the light from the nucleus 
and the total light of the object, setting the scene for the class 
of Seyfert galaxies that in summary show, bright, star-like nu- 


* Data displayed in Figs. [2] and 0 (Table [7] is an example) are only 
available in electronic form at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr 


cleus, and broad emissio n lines that cove r a w ide range of ion¬ 
isation states. iKhachikian & Weedmanl (l!974i) much later clas¬ 
sified Seyfert galaxies into types 1 and 2 depending on whether 
the HyS lines (Balmer series) are broader than the forbidden ones, 
or of approximately the same width. Typically, the width of for¬ 
bidden lines in both classes are of the order of “narrow emis¬ 
sion lines”, that is, of about 300-800 km s _1 , while the width of 
“broad emission lines” - including the hydrogen Balmer lines - 
is a bout 1000-6000 km s ~* ( Osterbrocklf19841) . 

iDavidson & Kinmaiil d 19781) found that MRK 359 appeared, 
to be lying at the low end of the line-width distribution and 
to have H/I and forbidden lines of comparable widths that 
were similar to those in Seyfert 2 galaxies (~ 300 km s _I ); 
furthermore, MRK 359 showed strong featureless continuum, 
and strong high-ionisation lines (e.g. [Fevn] and [Fex]); these 
properties are common in Seyfert 1 galaxies, which led to 
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the definition of yet another sub-class of AGN: the “narrow - 
line Seyfert 1” - hereafter NLSyl - galaxies dGaskelll 


lOsterbrock & Poggall985t lOsterbrock & PahariL 119831) . iKoskil 


( 1 978h and iPhillinsI (1 1 978h noted th a t MRK 42 showed_ 

ilar behaviour. lOsterbrock & Poggel (119851) studied eight such 
sources and concluded that they are characterised by (a) unusu¬ 
ally narrow Hfj lines, (b) strong Fen emission, (c) normal lu¬ 
minosities, and fd) H/i slightly weaker than in typical Seyfert 
1 galaxies. Conventionally, today sources are categorised as 
NLSyl if they show (a) a narrow width of the broad Balmer 
emission line with a FWHM(H/?) < 2000 km s_\ and (b) weak 
forbidden lines with rOml/15007/Hg < 3 (lOsterbrock & Pogg3. 

1 1 9851:iGoodrichl 11984Izhon et al 1 12006b . 

The first attempts to inve stigat e the radio prope rties of these 
systems were undertaken by lUlvestad et al.l d 1995b . among oth¬ 
ers, who studied seven NLSyls with the Very Long Array (VLA) 
in A configuration at 5 GHz. They found that the radio power 
was moderate (10 20 23 W Hz -1 ), the emission is compact (< 
300 pc), and in the two of three cases where radio axes could 
be found and high optical polarisation was detected, the radio 
axes were oriented perpendicularly to the electric vector posi¬ 
tion Angle (EVPA); the third one showed EVPA nearly parallel 
to the radio axis. Later. iMoranl d2000h also studied 24 NLSyls 
with the VLA in A configuration to confirm that most of the 
sources were unresolved and that they showed relatively steep 
spectra. fStenanian et al.l (12003 ) investigated 26 NLSyl galaxies 
and found that onl y 9 were detected in the FI RST catalogue 
( White et all 1 1 9971) . and all were radio-quiet. iKomossa et alJ 
(12006}) studied 128 NLSyls in a dedicated search for radio- 
loud (RL) NLSyls and concluded, among others, that (a) their 
morphology is sim ilar to compact steep-spectrum sources (CSS; 
iGallo et all I2006L discussed PKS 2004-447 as an example of 
NLSyl that was also CSS); (b) the radio-loudness R - defined 
as the ratio of the 6 cm flux to the optical flux at 4400 A 
(iKellermann et al.LI 1989b - is distributed smoothly up to the crit¬ 
ical value of 10 and covers about four orders of magnitude; (c) 
almost 7 % of the NLSyl galaxies are formally RL, but only 
2.5 % of them exceed a radio index/? > 100; (d) most RL NLSyl 
are compact steep-spectrum sources, than accreting close to or 
above the Eddington luminosity, Le ddi (e) the black-hole masses 
are generally at the upper observed end for NLSyl, but are still 
smaller than what is typically seen in other RL AGNs. 

Even before the detection of the RL subset of this class at 
high energies, it was clear that the y comprise a special source 
type for a number of reasons. lYuan et ahl d2008l) studied a sam¬ 
ple of 23 NLSyls with a radio-loudness exceeding 100 to find 
that these RL AGNs may be powered by black holes of mod¬ 
erate masses (10 6 -10 7 Mq) accreting at high rates and that 
they show a variety of radio properties reminiscent of blazars. 
IKomossa et al.l(l2006b pointed out that NLSyls provide an excel¬ 
lent probe for studying the physics scaling towards lower black- 
hole masses given their systematically low black-hole masses in 
the range 10 6 —10 s M e . Despite the fundamental differences be¬ 
tween RL NLSyls and blazars - in terms of black-hole masses 
and accretion disc luminosities (0.2-0.9 LEdd) - their jets seem 
to be behavingsimilarly and share the same properties (see also 
lFoschiniLl2012bl) . That alone has implications on the energy pro¬ 
duction and dissipation at_different scale s. 

According to IKomossa et al.l (i2006l) . an important question 
arising from the behaviour of these systems is whether or not the 
black-hole mass and radio-loudness are related and especially 
whether there is a limiting black-hole mass above which objects 
are preferentially RL, and whether or not RL and radio-quiet 
galaxies show the same spread in their black-hole masses. On 


this ground. iKomossal (' 2008 1 argued that the study of NLSyls 
needs to (a) include the radio and infrared properties of NLSyl 
galaxies when correlation analyses are performed, (b) determine 
the sufficient and necessary conditions for the onset of NLSyl 
activity, and (c) the investigation of whether the low black-hole 
mass is enough to explain the typically observed NLSyl charac¬ 
teristics. 

The first report for the detection of y-ray emission from 
a sourc e clas sified as an RL NLSy l was given bv lAbdo et al.l 
(l2009al) and iFoschini et al.l (1201 Oil , who discussed the detec¬ 
tion of significant GeV emission by the Fermi /LAT instrument 
from PMN J0948+0022. Already in its firs t year of operations, 
Fermi/ LAT detected a total of four NLSyls (lAbdo et all 120093) . 
As of today, there are seven RL NLSyls detected in the MeV- 
GeV ene rgy bands five of wh i ch with high significance (TS> 25, 
see also ID’ Ammando et ahl l2013bb IFoschini et all 12014b. A 
thoroug h review of the recent discoveries is given bv iFoschirvil 

(12012bl) . 


The multi-w avelength c a mpaign s that followed these 


discoveries (e . g. Abdoetal. . 2009b [at Gir oletti et all 1201 ll : 


IFoschini et all 2012 : Fuhrmann et al. 712 01 lb a nd the study of 


subsequently discovered NLSyls ( Abdo et all 12009(3) showed 
a clear blazar-like behaviour, indicating the existence of a rela¬ 
tivistic jet viewed at small angles. 

Despite the research activity that their high-energy detection 
has motivated, several questions remain open. Especially the ex¬ 
act nature and properties of NLSyls and the properties of the jet 
that seems to be present. Particularly in the radio regime, where 
a jet is expected to be dominant, the studies conducted so far 
have mostly been based on non-simultaneous datasets. 

In the current work we quantify some of the properties of 
the assumed radio jet emission. The motivation for this is to 
compare them with those seen in other classes of AGNs - es¬ 
pecially blazars. To do this, we studied the radio behaviour of 
four RL NLSyls detected in y rays by Fenni/LAT through regu¬ 
lar multi-band, single-dish radio monitoring with the Effelsberg 
100 m and IRAM 30 m telescopes. We focus on the variability 
properties - especially brightness temperatures and variability 
Doppler factors, which will (in a later publication) be combined 
with long-baseline radio interferometric measurements to esti¬ 
mate the viewing angles and Lorentz factors (Fuhrmann et al., 
in prep.). The energetics of the observed radio outbursts and the 
variability patterns of the broadband radio spectra are seen as in¬ 
dicators of the variability and emission mechanism at play. Very 
early re sults have been presented in several other studie s, in¬ 
cluding D^AmmandoeTal] (l2012b . lFuhrmann et al.l (1201 lb . and 
IFoschini et al.l ( 2012 ). 

Here we present the longest-term, multi-frequency radio 
monitoring datasets of the known y-ray-detected and radio- 
loud NLSyls to date. As an example, the 4.85 GHz light-curve 
lengths, range between longer than 1.5 years and longer than 5 
years. The monitoring was conducted at eight different frequen¬ 
cies for the two brightest sources and at six frequencies for the 
faintest. 

The paper is structured as follows: in Sect. [2] we review some 
archival information about the sources in our sample that is nec¬ 
essary for the following discussion. After describing the observ¬ 
ing methods in Sect. [3] we continue in Sect. [4] with the phe¬ 
nomenological description of the obtained light curves. The light 
curves discussed there comprise the basis for the flare decompo¬ 
sition method presented in Sect. [5] Subsequently, we continue in 
the frequency domain by studying the radio spectra in Sect. [6] 
In Sect. |7| we compute the radio powers for the 14.60 GHz 
datasets, while in Sect. [8] we discuss our polarisation measure- 
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ments. Summarising comments are discussed in Sect. [8] where 
we describe the general connection of all our findings. The main 
findings of our study are presented in Sect. [TO] which concludes 
the paper. 

Throughout the paper we assume a AC DM cosmology w ith 
H 0 = 70 km s -1 Mpc~' and = 0.73 ( Komatsu et al .. l201 ll) . 


2. The sample: what is already known 

We here investigate a sample of four sources, which are 

1. classified as NLSyls, 

2. radio-loud, 

3. detected at y rays by the FermifLAI, and 

4. satisfy certain observational constrains such as declination 
limits and a sufficiently high flux density to enable quality 
measurements at the Effelsberg (> 0.1 Jy) and IRAM tele¬ 
scopes (> 0.3 Jy). 

We briefly summarise what is already known of these sources. 
Table[T|summarises their positions, redshifts, and classifications. 

2. 1. J0324+3410 (1H 0323+342) 

The source J0324+3410 is of special interest mostly because of 
the morphology of its host galaxy and the estimated black-hole 
mass. From two HST archival images (each of 200 s exposure), 
IZhou et akl ( 2007 ) concluded that it is hosted by a one-armed 
spiral galaxy. The authors showed that its radio and X-ray emis¬ 
sion can be explained by jet synchrotron radiation, while the in¬ 
frared and optical light is dominated by thermal emission from a 
Seyfert nucleus. lAnton etliil d2008l) conducted B- and /(-band 
observations with the Nordic Optical Telescope (NOT). They 
suggested that its host resembles the mo rphology seen in the in¬ 
ner parts of Arp 10 found by ICharmandaris & Appletoni (119961) . 
implying that J0324+341 0 may b e a mer ger remnant. Following 
the method discussed by iGreene & Hoi d2005l) , they estimated 
the central black-hole mass to be 10 7 Mq; this value lies in the 
overlapping region between the black-hole mass distributions for 
N LSyls and bla zars and is similar to the value published earlier 
by IZhou et al.l d2007l) . 

Although the only case of an RL, Fermi- detected NLSyl 
for which the host is well resolved, J0324+3410, challenges the 
assumptions that powe rful relativistic jets form only in giant 
ellipti cal galaxies (e.g. iBottcher & Derrneil l2002t ISikora et all 
l2007h. 

Abdo et al. (2009c ) presented a thorough study including 
spectral energy distribution (hereafter SED) fits. The presence 
of a jet was already apparent in these studies; but in an episodic 
manner and not as a fully developed constantly broadband- 
emitting jet. The computed jet power placed the source in the 
BL Lac range (see jet power computations in Sect. [TJ. It is worth 
noting that they computed accretion rates that reach extreme val¬ 
ues of up to 90 % of the Eddington luminosity, yet another pe¬ 
culiar property of the RL, y-ray-loud NLSyls. 

2.2. J0849+5108 (SBS 0846+513) 

J0849+5108 was identified as an NLSyl galaxy by IZhou et al.l 
( 2005 ). They found clear evidence for emission from a rela¬ 
tivistic jet and a stellar component. Additionally, the emission 
lines show characteristics that classify it as a typical NLSyl 
(FWHM(H/3)^ 1710 km s _1 , [O iiiR, 0 07 =* 0.32 H J3 and strong 


Fen emission). The fact that_the source was previously classi¬ 
fied as a BL Lac object (Arp et al.. 1979 ) makes the case es¬ 
pecially interesting, possibly providing a handle on the link be¬ 
tween B L Lac objects and NLSyls. 

lYuan et all ( 2008 i) included the source in a very thorough 
study of a sample of 23 RL NLSyl galaxies and reported a black- 
hole mass of about 10 7 ' 4 Mq. On the basis of the detected signif¬ 
icant optical polarisation, they claimed that it has a jet. Finally, 
they found that during the high-energy states the optical con¬ 
tinuum is featureless, while at low states strong emission lines 
become obvious; this suggests that the source maybe a transition 
between a quasar and a BL Lac that displays its latter character 
only at flaring states. 

IFoschinil J2011bh was the firs t to re port its detection at 
y rays. Later, D’Ammando et al.l (1201 2l) discussed a flaring 
episode around June-July 2011. The multi-frequency datasets 
collected dur ing the campaigns following the detection are partly 
discussed in ID’ Ammando et all d2013cl) . Us ing the Very Long 
Baseline Array (VLBA) at 5, 8.4 and 15 GHz lD’Ammando et al.l 
( 2012 ) resolved the otherwise compact source to reveal a core¬ 
jet structure. The feature attributed to the jet shows a steep spec¬ 
trum. The core has been decomposed into two compact com¬ 
ponents. They discussed the detection of an apparent speed of 
about 8 c, indicating a relativistic jet. In summary, the power out¬ 
put (isotropic y-ray luminosity) of 10 48 erg s _1 on a daily scale 
(similar to that of luminous flat-spectrum radio quasars, FSRQs), 
the apparent superluminal motions and the radio variability ac¬ 
companied by spectral evolution indicate blazar-like relativistic 
jet._ 

ID’ Ammando et alJ d2012l) also considered the suggesti on of 
iGhisellini & Tavecchiol d2o'()8h and iGhisellini et alJ (120 1 ll) that 
the transition between FSRQs and BL Lac can be interpreted in 
terms of different accretion rates. They found that the position 
of the source in a y-ray photon index against luminosity plot 
(r y - L y ) also places it in the typical blazar territory. 


2.3. J0948+0022 (PMNJ0948+0022) 

J0948+0022 shows the usu al ch aract eristics of an NLSyl op¬ 
tical spectrum dWilliams et all 12002 ) with a FWHM(H/J) of 
1500±55 k m s" 1 , O m/H/5 < 3, and strong optical F n emission 
dZhou et all 120031) Although NLSyls are usu ally radio quiet 
(e.g. lUlvestad’etan [T995t iKomossa et all 120061) . J0948+0022 
was found to be th e first very ra dio-loud NLSyl wi th a R > 10 3 
( Zhou et al ., 120031) . In addition. IZhou et al. (20031) reported an 
inverted radio spectrum and high brightness temperatures, which 
strongly supports that there is a relativistic, Doppler-boosted jet 
seen at small viewing angles. The source has also been the first 
NLSyl detected at high-energy y rays by Fermi /LAT during its 
first months of operation. Fermi/LAI confirmed the relativistic 
jet emission from this source and established NLSyls as a new 
class of y-ray emitting AGNs dFoschini et al.Ll2010tlAbdo et al .. 
l2009allbl) . 

At pc scales, the source appears as one-sided VLBI structure 
dominated by a compact (< 55 r tas), b right central component 
dDoi et al.L l2006t iGiroletti et all 2 01 lb. Kpc -scale radio emis¬ 
sion has also been found bv iDoi et al.l d2012l) with a two-sided 
extension of the core and a northern extent of 52 kpc in pro¬ 
jected distance. From these studies a relativistic all y boos ted jet 
with Doppler factors D > 1 (e.g. D > 2.7 - 5.5 Doi et all 12006 ) 
has been inferred, in agreement with Doppler boos ting inferred 
from previous variability studies (e.g. iFuhrmann et all 1201 ll) . 

This NLSyl is variable across the whole electromagnetic 
spectrum. Intense multi-wavelength campaigns performed over 
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Table 1. The four monitored NLSyls with their observed positions, redshift, and classification. 


Source ID 

Survey name 

RA 

(hh:mm:ss.s) 

DEC 

(dd:mm:ss.s) 

z 

Class 

J0324+3410 

B20321+33B 

03:24:41.2 

+34:10:45.1 

0.0610“ 

NLSyl' 





0.0629 b 


J0849+5108 

SBS 0846+513 

08:49:58.0 

+51:08:29.0 

0.584701 c 

NLSyl 2 

J0948+0022 

PMN J0948+0022 

09:48:57.3 

+00:22:25.6 

0.585102 c 

NLSyl 3 

J1505+0326 

PKS 1502+036 

15:05:06.5 

+03:26:31.0 

0.407882 c 

NLSyl 4 


References, (a) IMarcha et all 119961): ( b) IZhou et all 120071) : (c) iMarcha et al.l (Il99fih : (1) IZhou et all 120071) : (2) iGreene & Hoi d2005h : (3) 
EE ou et a 0 120031) : 141 IZhou et al.ll200fil) 


the past years revealed continuou s act ivity with an exceptional 
outburst in 2010 dFoschini et al.L I20TT1 ) that occurred at all spec¬ 
tral bands. Together with detailed SED modelling and the detec¬ 
tion of polarised emission, these findings confirm the powerful 
relativistic jet in the source seen at small viewing angles. 


2.4. J1505+0326 (PKS 1502+036) 

Based on its optical characteristics (FWHM(H/5) = 1082 ± 
113 km s -1 , [O iii|/H/j~ 1.1, and a strong o ptica l Fe n bum p), 
this source was classified as NLSyl (e.g. lYuan et al.L 120081) . It 
exhibits one of the highest radio-loudness parameters among 
NLSyls (R = 1549, see also Sect. [ 53 ). Early VLBI ob¬ 
servations revealed a compact source (marginally resolved at 
milli -arcseco nd scales ) with an inverted radio spectral index 
(iDallacasa et al.Lll998i) . 

lAbdo et ak! (i2009d) reported the first detection of y -ray emis¬ 
sion from the source by FermifLAT. From SED modelling the 
authors inferred jet powers similar to those in J0948+0022 and 
those typically observed in powerful quasars. A detailed multi¬ 
wavelength campaign carried out between 2008 and 2012 did not 
reveal any significant variability at y rays, in contrast to promi- 
nent flux density and spectr al variability seen in the radio regime 
dD’ Ammando et all . 12013 ah . The 15 GHz VLBI imaging showed 
a one-sided core-jet structure on pc scales. No significant proper 
motion of jet components was detected. The radio variability and 
VLBI findings together with an inferred high apparent isotropic 
y-ray luminosity strongly support a relativistic, Doppler-boosted 
jet in this case as well. 


2.5. Radio-loudness 


Of the four sources in our study, J0948+0022, J0849 +51 08, and 
J1505+0326 are very RL, with/^i 4 = 194-793 llKomossa et all 
120061) . 1445 and 1549 (IYuan et all l2008h. res pective ly; where 


R is the radio index according to Kellermann et al. (19891). 


J0324+3410 is only mildly RL with R$ = 38 - 71 (izhou et all 
[2007). We note that each value of R is uncertain by a factor of a 
few to 10 , reflecting uncertainties in extinction correction, opti¬ 
cal host contribution, and variability in the radio or optical bands. 
The monitoring presented here, reveals variability in the radio 
bands of up to several ten of percent (see Table [ 8 ), but not high 
enough to modify the previous source classification as RL or 
radio-quiet. The values of the radio indices reported in the liter¬ 
ature are therefore adequate for the scope of this paper and are 
shown in Table [2] 


Table 2. Radio-loudness of the sources in our sample. The mod¬ 
erate variability cannot cause remarkable changes in the R in¬ 
dices, making the published values representative enough. 


Source 

Radio Index (a) 

Reference 

J0324+3410 

R s = 38 - 71 

(1) based on HST optical flux 

J0948+0022 

m 

ON 

r- 

1 

ON 

II 

(2) based on USNO and GSC 
catalogues (note the variable 
optical flux) 

J0849+5108 

R l 4 = 1445 

(3) based on SDSS optical flux 

J1505+0326 

Ria = 1549 

(3)... 


Notes. With R 5 = S 5 ghz/Sb. Ria = S 1 .4ghz/5'44ooa Kellermann’s 
radio index. Under the assumptions of Kellermann et al.: R\ .4 = I. 9 R 5 . 


References. (1) 


f3) 1 Yuan et all ( 2008) 


Zhou et all 120071) : (2) [Komossa et afl 120061) : 


3. Observations and data reduction 

The light curves and radio SEDs presented here have been 
observed in the framework of the F - GAM MA programme 
dFuhrmann et all 2001 : Angelakis et all 12010 ). They cover the 
frequency range from 2.64 to 142.33 GHz. Below 43.05 GHz 
the measurements were conducted with eight different hetero¬ 
dyne receivers mounted on the secondary focus of the 100 m 
Effelsberg telescope. The observations at 86.24 and 142.33 GHz 
were obtained with the 30 m IRAM telescope. In Table [4] their 
most important operational characteristics are summarised . 


3.1. Effelsberg observations 

The Effelsberg station covers the band between 2.64 and 
43.05 GHz. The systems at 4.85,10.45 and 32 GHz are equipped 
with multiple feeds allowing differential measurements meant to 
remove atmospheric effects (e.g. emission or absorption fluctua¬ 
tions). The rest are equipped with only single feeds. All receivers 
have circularly polarised feeds. The observations are made in 
cross-scanning mode, that is monitoring the telescope response 
while driving over the source position in two different orthog¬ 
onal directions (in our case, azimuth and elevation). Necessary 
post-measurement corrections include the following: 

1. Pointing-offset correction: meant to correct for the power 
loss caused by possible divergence of the actual source posi¬ 
tion from that observed by the telescope. This effect is of the 
order of a few percent. 

2. Atmospheric-opacity correction: meant to correct for the 
attenuation caused by the transmission through the terres- 
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Table 3. Flux densities of the standard calibrators used at 
Effelsberg. 


Source: 

3C48 

3C161 

3C 286 

3C295 

NGC 7027" 

S 2.64 

9.51 

11.35 

10.69 

12.46 

3.75 

S 4.85 

5.48 

6.62 

7.48 

6.56 

5.48 

S 8.35 

3.25 

3.88 

5.22 

3.47 

5.92 

S 10.45 

2.60 

3.06 

4.45 

2.62 

5.92 

S 14.60 

1.85 

2.12 

3.47 

1.69 

5.85 

S 23.05 

1.14 

1.25 

2.40 

0.89 

5.65 

S 32.00 

0.80 

0.83 

1.82 

0.55 

5.49 

S 43.00 

0.57 

0.57 

1.40 

0.35 

5.34 


Notes. (a) The flux density of NGC7027 is corrected for beam- 
extension at frequencies above 10.45 GHz. 


References. Th e flux densi t ies of the calibrato r s are taken from 
lOtt et all dl994l) . iBaars et all dl977l) . IZiilstra et all d2008f) . and Kraus 
priv. comm. 

trial atmosphere. This effect can be significant especially at 
higher frequencies where the atmospheric opacity becomes 
significantly high. 

3. Elevation-dependent gain correction: correcting for sensitiv¬ 
ity losses caused by small-scale departures of the primary 
reflector’s geometry from that of an ideal paraboloid (as ex¬ 
pected for a homology-designed reflector). The magnitude 
of this effect is constrained to within a few percent. 

4. Absolute calibration: converting the measured antenna tem¬ 
peratures to SI units by reference to standard candles, that 
is, calibrators. The standard candles used for our programme 
along with the flux densities assumed for them are shown 
in Table [3] Specifically, in the case of NGC7027, a partial 
power loss caused by its extended structure relative to the 
beam size above 10.45 GHz was corrected for. 

The five-year mean magnitudes of these effects as they are 
observed by the F-GAMMA programme are summarised in 
Tabled 

3.2. IRAM observations 

Observations with the IRAM 30 m telescope were made 
within the F-GAMMA monitoring programme and the more 
general flux monitoring conducted by IRAM (Institut de 
Radioastronomie Millimetrique, lUngerechts et all Il998l) . Data 
from both programmes are included in this paper. The observa¬ 
tions were conducted with the newly installed EMIR receiver 
( Carter et all 12012 ) using the 3 and 2 mm bands (each with lin¬ 
ear polarisation feeds) tuned to 86.24 and 142.33 GHz and the 
narrow-band continuum backends (1 GHz bandwidth) attached. 
Observations of J0324+3410 and J0948+0022 were performed 
with cross-scans in azimuth and elevation direction and wobbler¬ 
switching along azimuth with a frequency near 2 Hz. Each cross¬ 
scan was preceded by a c alibration scan to obtain in stantaneous 
opacity information (e.g. lMauersberger et all 11989 ). After data 
quality control, the sub-scans of each scanning direction were 
averaged and fitted with Gaussian curves. In the next step, each 
amplitude was corrected for small remaining pointing offsets and 
systematic gain-elevation effects. This operation has an effect of 
about 1 % at 86.24 GHz and 4 % at 142.33 GHz (mean point¬ 
ing offsets are 1.8 "for both receivers). The latter correction, and 
given the IRAM 30 m elevation-dependent gain curves, amounts 
to < 5 %. The conversion to the standard flux density scale was 


Table 5. The average fractional effect of each post-measurement 
correction applied to the data for each observing frequency. The 
numbers represent the 5-year average behaviour as it is observed 
by the F-GAMMA programme. 


Frequency 

(GHz) 

Pointing 

Correction 

(%) 

Opacity 

Correction 

(%) 

Gain 

Correction 

(%) 

2.64 

Effelsberg 

0.5 2.3 

0.0 

4.85 

0.4 

2.5 

1.2 

8.35 

0.5 

2.5 

0.9 

10.45 

1.2 

3.1 

1.2 

14.60 

1.3 

2.9 

1.3 

23.05 

1.6 

8.1 

2.1 

32.00 

3.1 

7.9 

2.9 

43.05 

5.1 

20 . 

2.1 

86.24 

IRAM 

1.0 

5.0 

142.33 

4.0 


5.0 


made using frequent observations of primary (Mars, Uranus) and 
secondary (W30H, K3-50A, NGC7027) calibrators. 

3.3. Error estimates 

In every operational step the errors were propagated formally 
assuming Gaussianity. The exact details are discussed in sepa¬ 
rate papers (Nestoras et al. in prep., for the IRAM observations; 
Angelakis et al. in prep., for the Effelsberg observations). An 
indicative empirical measure of the uncertainty in a measure¬ 
ment can be given by the fractional fluctuations seen in sources 
known to be intrinsically stable, that is to say, the calibrators. 
The datasets available for these sources cover the entire period 
of observations and include - cumulatively - all possible sources 
of fluctuations. This variability - seen in sources intrinsically 
stable - must also be present in the light curves of the targets. 
In Table[6]we quote the mean flux density at each frequency for 
each calibrator used at Effelsberg and NGC 7027 used at IRAM. 
Note that the flux density for each such source was computed 
using the mean calibration factor of each session. There we also 
show the modulation index defined as m — 100 ■ with cr being 
the standard deviation and (S ) the mean flux density, m remains 
at levels of a few percent even at frequencies above 23.05 GHz 
where the troposphere becomes very disturbing. 

3.4. Cross-telescope calibration 

Because below we discuss the dynamics of radio SEDs observed 
partly with the Effelsberg and partly with the IRAM 30 m tele¬ 
scope, it is essential to address the cross-telescope calibration 
accuracy that might introduce artefacts in the observed spectral 
shapes and their apparent temporal behaviour. An empirical yet 
reliable measure of its goodness can be provided by NGC 7027, 
which exhibits a flux density high enough to be detected by both 
instruments with high signal-to-noise ratios (S/N) and repeat¬ 
edly. NGC 7027 has a very well defined and analytically de¬ 
scribed time-dependent convex spectrum (IZiilstra et al.L 120081) 
that at frequencies above roughly ~ 10 GHz can be approxi¬ 
mated by a power law of the form S oc y 01 . In Fig. Q] we show 
the flux densities averaged over the entire period discussed here 
and for all frequencies. Each circle denotes an average flux den- 
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Table 4. Receiver characteristics. 


V 

(GHz) 

A 

(mm) 

Av 

(GHz) 

T S y S 

(K) 

FWHM 

(") 

Feeds 

Polarisation 

Aperture Efficiency 
(%) 

2.64 

110 

0.1 

17 

260 

1 

LCP, RCP 

58 

4.85 

60 

0.5 

27 

146 

2 

LCP, RCP 

53 

8.35 

36 

1.2 

22 

82 

1 

LCP, RCP 

45 

10.45 

28 

0.3 

52 

68 

4 

LCP, RCP 

47 

14.60 

20 

2.0 

50 

50 

1 

LCP, RCP 

43 

23.05 

13 

2.7 

77 

36 

1 


30 

32.00 

9 

4.0 

64 

25 

7 

LCP 

32 

43.05 

7 

2.8 

120 

20 

1 


19 

86.24 

3 

8.0(1 GHZ used) 

-65“ 

29 

1 

HLP, VLP 

63 

142.33 

2 

4.0(1 GHZ used) 

-65“ 

16 

1 

HLP, VLP 

57 


Notes. The entries in each column is as follows: 1: central frequency, 2: receiver bandwidth, 3: system temperature, 4: sensitivity, 5: full width 
at half maximum ( FWHM ), 6: number of available feeds, 7: available polarisation channels, 8: telescope effective area at the corresponding 
frequency. 

(|I) The values quoted for the IRAM receivers are typical values for the receiver temperatures T IX and not system temperatures, hence they do not 
include atmospheric contributions, background emission, etc. 


Table 6. The mean flux density and modulation index of the calibrators at each frequency. 


Source 

Observable 

Units 

2.64 

4.85 

8.35 

10.45 

14.60 

23.05 

32.00 

43.05 

86.24 

142.33 

3C286 

N 


102 

108 

108 

111 

106 

101 

89 

46 




(S') 

(Jy) 

10.710 

7.477 

5.2189 

4.453 

3.476 

2.408 

1.840 

1.429 




m 

(%) 

0.9 

0.5 

0.7 

1.1 

1.9 

2.9 

3.5 

4.4 



3C48 

N 


88 

91 

92 

92 

93 

83 

72 

28 




<s > 

(Jy) 

9.579 

5.509 

3.267 

2.613 

1.872 

1.162 

0.817 

0.599 




m 

(%) 

0.7 

0.6 

0.9 

1.2 

1.7 

2.9 

4.6 

4.5 



3C161 

N 


30 

35 

29 

30 

30 

22 

19 

4 




<s> 

(Jy) 

11.361 

6.597 

3.793 

2.961 

2.053 

1.221 

0.822 

0.597 




m 

(%) 

1.0 

0.8 

1.4 

2.2 

3.3 

5.4 

7.3 

2.3 



3C295 

N 


38 

38 

34 

35 

35 

25 

24 





<s> 

(Jy) 

12.499 

6.546 

3.434 

2.588 

1.690 

0.880 

0.557 





m 

(%) 

0.9 

0.8 

1.3 

1.5 

2.6 

6.4 

6.7 




NGC7027 

N 


64 

73 

79 

75 

71 

62 

52 

25 

61 

54 


(S) 

(Jy) 

3.693 

5.449 

5.903 

5.912 

5.738 

5.455 

5.243 

5.146 

4.920 

4.611 


m 

(%) 

1.2 

0.5 

1.1 

1.2 

2.1 

3.8 

5.3 

5.3 

2.7 

4.6 


sity. The red symbols (with v > 10.45 GHz) are the Effelsberg 
measurements that were used for the model fit, and the red line 
is the result of the fit. The grey area is confined by 1 cr of each av¬ 
erage flux density. Extrapolating the fitted spectrum (red dotted- 
dashed line) towards the higher IRAM band and comparing these 
values with the measured IRAM 30 m flux densities yields dif¬ 
ferences better than 3 %, specifically, 2.7 % at 86.22 GHz and 
2.8 % at 142.33 GHz), indicating a high-quality cross-telescope 
calibration. 


3.5. Analysis methods 


The current section introduces notions and methods that are used 
below, even if they are repeated briefly in the corresponding sec¬ 
tions. 


Internal shocks causing variability: One of the aims here 
is to study the mechanism that may be causing the observed 
variability. Throughout the text it is assumed that this could 
well be caused by internal shocks that propagate in the jet, im¬ 
printinga_s£ecific_signaturein the_£adio_Jightcurves. In this 
model dMarscher & GeaA 1 1 9851 iTtirler et~akl l2000l) . the syn¬ 
chrotron self-absorbed component is undergoing distinct evo¬ 
lutionary stages, each characterised by a different energy-loss 
mechanism. The followed path then imprints a distinct phe¬ 
nomenology on the radio SEDs, making this model easily quan¬ 
tifiable and testable. It has been argued that the implementation 
of this scenario in a system with simply a steep-spectrum qui¬ 
escent jet is en ough to reproduce the o bserved plurality of phe- 
nomenologies (lAngelakis et all 12012a) . 


The variability brightness temperature: As discussed exten¬ 
sively in Sect. 0 the variability brightness temperature is a mea- 
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Fig. 1. Observed radio SED of NGC 7027 dZiilstra et all 120081) 
over the same period of time as the time baseline covered here. 
All the Effelsberg frequencies from 2.64 up to the highest IRAM 
frequency of 142.33 GHz are shown. The filled circles (red 
and black) denote the measurements. The red symbols mark 
the Effelsberg measurements that where used in the fit, the red 
dotted-dashed line is the fitted spectrum. The grey area denotes 
the lcr region around each data point. The agreement of the ex¬ 
trapolated values with those measured by IRAM is better than 
3 %. 


Table 7. Example of a light-curve file as published online at 
CDS. The current extract is for J0324+3410 at 4.85 GHz. 


MJD 

5 

(Jy) 

Error 

(Jy) 

55408.300 

0.485 

0.006 

55458.148 

0.395 

0.007 

55487.059 

0.440 

0.011 

55515.026 

0.423 

0.005 

55570.904 

0.374 

0.005 


in TableflOl In Sect. [5] we explain that the limiting value i s bas ed 
on the equipartition argument as proposed bv lReadheadl (1 1994b . 

4. Radio light curves 

The dataset discussed here is the result of the longest-term multi¬ 
frequency radio monitoring for RL and y-ray loud NLSyls avail¬ 
able. It covers the period between early-2009 and mid-2014. In 
Figs.[2]and[3]we present the available light curves at all observ¬ 
ing frequencies. We only show measurements with a signal-to- 
noise ratio better than 3. All the light curves are also available 
online at the CDS. Table[7]is an example of an online light-curve 
file for J0324+3410 at 4.85 GHz. It contains the following infor¬ 
mation: Column 1 lists the MJD of a measurement. Col. 2 gives 
the flux density in units of Jy, and Col. 3 lists the associated error. 


sure of the energetics of the associated event. It is computed on 
the basis of the light travel-time argument and depends on the 
magnitude of the flux density variation, SS and the time span 
needed for that variation, St. The flare decomposition method 
(Sect. [5} aims at separately estimating exactly those parame¬ 
ters for each event. The variability brightness temperature at the 
source rest frame in K, is given by 


T vai = 1.64- 10 10 


SS ■ dl ■ A 2 
St 2 ■( 1 + z) 4 ’ 


(1) 


where 


SS is the increase in flux density (at the observer’s 
frame) in units of ly, 

St is the time span needed for that increase (at the ob¬ 
server’s frame) in units of days, 
d L is the source luminosity distance in units of Mpc, 

A is the observing wavelength in units of cm, and 
z is the source redshift. 


all measured at the observer’s frame. 

Any excess from an independently computed intrinsic limit 
is then interpreted in terms of the Doppler-boosting factor D as 

Z) = (1 +z) 3 a/S’ (2) 

where 

a is the source spectral index 1 with Xkv", 

T le f is the limiting value of the brightness temperature. 


The mean sampling - averaged over all frequencies - ranges 
from one measurement every 24 days for J0324+3410 and 
J0948+0022 (approximately the cadence of the most commonly 
used calibrator, 3C286), and 29 days for J1505+0326, to more 
than 100 days for J0849+5108. Except for 10849+5108, for 
which the dataset does not cover a substantial number of ac¬ 
tivity cycles, the remaining three sources show intense variabil¬ 
ity at practic ally all frequencies, a behaviour commonly seen in 
blazars le.g. iRichards et all 1201 lL lAller et akl 1201 it iBoettcheil 
120 1 2t iBCttcherL 12004b . The phenomenologies seen in these light 
curves vary significantly and mostly in terms of the observed 
amplitude of the different events. In these terms, J1505+0326 
shows very weak outbursts. 10324+3410 shows at least two ma¬ 
jor events at the high-end of the bandpass at around MJD 56300 
and 56600, which disappear at lower frequencies. J0948+0022 
shows frequent recursive events of activity that occur consider¬ 
ably fast; their characteristic times are of the order of 40-50 days 
at the highest frequencies (see Table IB.2l >. Most of the observed 
events can be identified at all frequencies (apart from the very 
low end of the bandpass) and with a time delay between frequen¬ 
cies (e.g. see event around MJD 56500 for J0948+0022), which 
is indicative of intense spectral evo lution th at is also routinely 
seen in blazar light curves (e.g. Rani et akl 12013 ). The variabil¬ 
ity amplitude is frequency dependent, following the typical fash¬ 
ion: higher frequencies vary more. For J0324+3410for example, 
the amplitude of the most prominent event is more than twice as 
large at 32.00 GHz as at 14.60 and more than six times lesser 
at 10.45 GHz; yet another indication that the variability mech¬ 
anisms jieem_to be the same as those acting in typical blazars 
(e.g. IValtaoia et all Il992h . Table [8] reports some characteristic 
parameters for these datasets. Specifically, for each frequency 


The value for the spectral index in the previous equation is cho- 1 Here we assumed that the emission comes from a single blob. For 
sen to be the mean values in the corresponding sub-band as given a continuous jet the index 3 + a should be replaced by 2 + a. 
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Fig. 2. Radio light curves available for J0324+3410 (left column) and J0849+5108 (right column) at all available frequencies. From 
top to bottom we present the light curves three different frequency bands: low: 2.64, 4.85 and 8.35 GHz, intermediate: 10.45, 14.60 
and 23.05 GHz, and high: 32, 43.05, 86.24 (when available) and 142.33 GHz (when available). At the very bottom the datasets are 
shown over-plotted together for comparison. For the same reason, for each source and axes the boundaries are kept identical. Lines 
connecting the data points have been used everywhere to facilitate visual inspection. Each frequency is consistently represented by 
the same colour and symbol. Only data points with a signal-to-noise ratio better than 3 have been used. 
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Fig. 3. Radio light curves available for J0948+0022 (left column) and J1505+0326 (right column) at all available frequencies similar 
as is Fig. [2] 
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and source we report the total number of data points N, the 
dataset length At, the mean flux density ( S ), the standard de¬ 
viation around that mean cr, and the modulation index m. 

In the following we roughly describe the most noteworthy 
phenomenological characteristics of the light curves separately 
for each source. 

4.1. J0324+3410 

The multi-frequency light curves of J0324+3410 are shown in 
Fig. [2] It shows a clearly different behaviour between high- and 
low-frequency bands despite the rather sparse sampling. Of the 
main events present in the 86.24 and 142.33 GHz light curves - 
around MJD 56040, 56300 and 56620 - the first one disappears 
already at 10.45 GHz, the second is still present at 8.35 GHz, and 
the last can be traced down to 4.85 GHz. These events also show 
different characteristics such as rise and decay times. Generally, 
the source activity is very moderate at the lowest frequencies. 

The neighbouring frequencies show associated events and 
strong evidence of intense spectral evolution. Differences be¬ 
tween similar frequencies can be regarded as unusual because 
the events disappear fast over frequency. At the highest frequen¬ 
cies, flux density variations by factors of about 5 or more are 
present, which are absent at the lowest frequencies. The spectral 
evolution may be an even better proxy of this acute behaviour, 
as we discuss later. 

Finally, we note that the baseline at all frequencies is prac¬ 
tically flat, an indication that all the variability incidents evolve 
fast and dissipate in a dominating relic jet that either does not 
display any signs of variability or does so at a very slow pace; 
the dynamics of the radio SED shape point toward such an inter¬ 
pretation, as well. 

4.2. J0849+5108 

The very few available SEDs for this source indicate that this 
is another interesting and very active source. Unfortunately, the 
lack of a large enough and adequately sampled dataset prevents 
any systematic quantitative analysis. Nevertheless, the source 
shows activity cycles at all frequencies,which evolve system¬ 
atically slower at lower frequencies (see the slow rising trend 
that becomes faster toward higher frequencies). It remains to be 
studied whether there is spectral evolution and what its char¬ 
acteristics are. Clearly, longer time-baselines are needed before 
any sensible results can be reached. 

4.3. J0948+0022 

J0948+0022 is among the best studied sources because it has 
been the first RL NLSyl to be detected by Fermi. The light 
curves shown in Fig. [3] clearly indicate intense and repetitive 
variability that is prominent at all frequencies except possibly 
at 2.64 GHz, where most of the outbursting events have already 
smeared out. The light curves show variability corresponding to 
factors of more than 3 even at frequencies as low as 10.45 GHz 
or below. As an example, the light curve at 8.35 GHz undergoes 
a series of clearly discernible outbursts at MJDs, approximately 
55020, 55240, 55490, 55725, 55985, and 56450. These events, 
which are single flares or sub-flares, appear in practically every 
densely enough sampled light curve except for that at 2.64 GHz 
(where the events have already disappeared). Even more interest¬ 
ingly, the moment of occurrence of an event appears at progres¬ 
sively later times as the frequency decreases, which indicates 


opacity effects. That is, a local maximum in a light curve cor¬ 
responds to the instant at which the emitting plasma radiation 
becomes optically thin at the observing frequency. If the flare 
is associated with the emergence of an adiabatically expanding 
plasmon, this instant should indeed appear at progressively later 
times as the frequency drops. This claimed spectral evolution is 
very clearly seen in Fig|4]as we discuss below. Finally, the pace 
at which the observed flux density ( dS/dt) increases during a 
flare is clearly a function of frequency, with higher frequencies 
showing much larger derivatives. 

4.4. J1505+0326 

The last source in our list of targets is J1505+0326. Its mean flux 
density of 507 mJy at 4.85 GHz, makes it the brightest member 
of the sample. Significant variability is also present, but with 
qualitatively different characteristics. The most important differ¬ 
ence - at least compared with J0948+0022 - is that the clear 
spectral evolution with the delay of the peak as a function of fre¬ 
quency is not as obvious here. Pairs of adjacent frequencies such 
as 4.85 and 8.35 GHz seem to show events in phase, although 
the lowest frequency, 2.64 GHz, does indeed lag by clearly dis¬ 
cernible time spans. At intermediate frequencies - for example 
at 10.45 GHz - the amplitude of variability reaches moderate 
factors of 1.5. Finally, frequencies below 10.45 GHz show a 
long-term decaying trend modulated by the faster variability dis¬ 
cussed above. 

5. Flare decomposition: parametrising variability 

The parametrisation of flux density outbursts - in terms of ampli¬ 
tude and time scale - is commonly used as a method to constrain 
the variability brightness temperature 7' var associated with the 
event, on the basis of causality arguments. The variability bright¬ 
ness temperature can only provide a lower limit of the intrinsic 
brightness temperature T b- Values of 7’ var in excess of indepen¬ 
dently calculated limiting values of the brightness temperature 
are generally attributed to Doppler boosting, which provides a 
handle for the computing limiting Doppler factors, D. Assuming 
an equipartition brightness temperature upper limit of 5 • 10 10 K 
( ReadheadLll99l . one can estimate the Doppler factors required 
to explain the observed excess. 

The combination of variability Doppler factors with Very 
Long Baseline Interferometry (VLBI) measurements of the ap¬ 
parent speeds allows compu ting the pla smo id bulk velocity and 
jet viewing angle (e.g. lLahteenmaki & ValtaoiaUl999l) . Here, we 
are equally interested in all these properties for our four RL 
NLSyls, and most importantly, in investigating the possible dif¬ 
ferences in the characteristics of different flares in the same light 
curve, rather than retrieving the characteristics of an average be¬ 
haviour. 

A practice commonly followed in variability studies is the 
implementation of time-series analysis methods that are de¬ 
signed to reve al su ch quantities; for example the structure func¬ 
tion anal ysis (| Simonetti et al.l.|l985h and the discrete correlation 
function dEdelson & Kroliki 119881) . One of the most important 
caveats of such methods, however, is that they are extremely sen¬ 
sitive to parameters that are difficult to determine in moderately 
sampled light curves such as the onset of a flaring episode or the 
shape of the temporal behaviour of the measured flux density. In 
fact, even minor changes in such parameters can result in differ¬ 
ences in the estimation of the variability brightness temperature 
beyond an order of magnitude. Furthermore, in most cases, these 
tools are designed to detect a dominant behaviour that smears 
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Table 8. Summary of the mean flux densities and corresponding standard deviations for each observing frequency and source 
presented in Figs. [2]and [3] For every entry we also report the number of data points N and the modulation index m = 100 ■ cr/ (S) 
as a measure of the apparent variability amplitude. Only data points meeting the condition of S/N> 3 are included. 


Source 

Observable 

Units 

2.64 

4.85 

8.35 

10.45 

14.60 

23.05 

32.00 

43.05 

86.24 

142.33 

J0324+3410 

N 


42 

45 

48 

48 

43 

14 

27 

8 

17 

17 


At 

(days) 

1310 

1310 

1310 

1310 

1310 

1260 

1260 

517 

899 

899 


(S) 

(Jy) 

0.461 

0.400 

0.380 

0.377 

0.378 

0.500 

0.432 

0.353 

0.541 

0.518 


cr 

(Jy) 

0.061 

0.064 

0.077 

0.083 

0.113 

0.230 

0.272 

0.115 

0.172 

0.164 


m 

(%) 

13 

16 

20 

22 

30 

46 

63 

32 

32 

32 

J0849+5108 

N 


11 

10 

12 

11 

9 

4 

8 

3 




At 

(days) 

981 

581 

1009 

581 

980 

552 

980 

152 




(S) 

(Jy) 

0.203 

0.214 

0.244 

0.266 

0.287 

0.363 

0.391 

0.529 




cr 

(Jy) 

0.012 

0.019 

0.028 

0.028 

0.040 

0.043 

0.103 

0.079 




m 

(%) 

6 

9 

12 

11 

14 

12 

26 

15 



J0948+0022 

N 


56 

63 

63 

65 

57 

31 

43 

14 

30 

23 


At 

(days) 

1863 

1863 

1863 

1863 

1862 

1739 

1862 

1766 

1644 

1502 


(S) 

(Jy) 

0.177 

0.255 

0.356 

0.390 

0.442 

0.555 

0.488 

0.569 

0.596 

0.530 


cr 

(Jy) 

0.039 

0.086 

0.142 

0.157 

0.173 

0.155 

0.188 

0.226 

0.193 

0.169 


m 

(%) 

22 

34 

40 

40 

39 

28 

38 

40 

32 

32 

J1505+0326 

N 


37 

42 

39 

43 

37 

15 

26 

7 




At 

(days) 

1204 

1261 

1261 

1261 

1261 

1233 

1261 

301 




<s> 

(Jy) 

0.456 

0.507 

0.503 

0.488 

0.465 

0.463 

0.408 

0.377 




cr 

(Jy) 

0.076 

0.085 

0.080 

0.072 

0.066 

0.050 

0.079 

0.103 




m 

(%) 

17 

17 

16 

15 

14 

11 

19 

27 




out possible significant differences in the characteristics of indi¬ 
vidual flares of the same source and even at different observing 
frequencies. 

To overcome such complications, we introduce a novel 
method for 

a. first creating “cumulative’' light curves by conveniently shift¬ 
ing and re-normalising the observed light curves. This oper¬ 
ation is meant to highlight the flares that are detectable at a 
wide range of frequencies, 

b. and subsequently subjecting those light curves to all neces¬ 
sary operations to extract the desired parameters (i.e. flare 
onset, duration, amplitude, etc.). 

The guiding principles while developing this approach were to 

1 . avoid complications introduced by the superposition of si¬ 
multaneously acting processes. For example, time-series 
analysis methods often return unrealistically long timescales 
only as the result of having - for example - a long-term 
almost-linear trend underlying much faster events. 

2 . to accommodate a generic approach in the treatment of ev¬ 
ery flare. That is, to parametrise each event independently 
(for each source and frequency) and investigate the possibil¬ 
ity of different behaviours (and possibly variability mech¬ 
anisms) acting in the same source at different times. For 
J0324+3410, for example, the most prominent event seems 
to demand such an approach because its phenomenology is 
very different from the rest (Fig. ED 

All the details of the method are discussed in AnnendixlAl 

5.1. Results of the flare decomposition 

Since the variability brightness temperature 7' vilI comprises only 
a lower limit of the brightness temperature 7 b, the higher the 7 var 


estimate, the better 7 b is constrained. Henceforth, from all the 
different estimates of 7 var we compute in the following analysis, 
the highest value is regarded as the most meaningful one. 

The method described above was applied to J0324+3410, 
J0948+0022, and J1505+0326, for which sufficiently long 
datasets were available and which showed a significant num¬ 
ber of activity “cycles” (5, 6 and 4, respectively). Table [9] sum¬ 
marises the results of this analysis. The highest frequencies are 
missing clearly because of the lack of data points at these bands. 

5.1.1. Flare decomposition for J0324+3410 

Table IB. 21 summarises the results of the flare decomposition 
method for all the significant events seen in the source light- 
curves shown in Fig. [2] As significant events are regarded flux 
variations that are above the noise level and can be detected 
at several wavelengths. As can be seen there, the variability of 
J0324+3410 is characterised by five fast variations shorter than 
66 days (even at the lowest frequencies) with a relatively low am¬ 
plitude; the largest amplitudes of the flares is generally reached 
at 14.60 GHz. 

The flare occurring at around MJD 56313 shows exceptional 
characteristics. Its maximum is reached at 32.00 GHz, and at that 
frequency its amplitude is around ten times larger than the aver¬ 
age of the other flares. The extraordinary phenomenology of this 
event alone could justify the introduction of an alternative anal¬ 
ysis method like the one described here. A classical Structure 
Function analysis, for example, would have smeared the event 
out. Below 10.45 GHz it moderates its behaviour and displays 
amplitudes comparable to that of other events. Its temporal be¬ 
haviour is similar to that of the immediately subsequent event 
which peaks around MJD 56356 at 32.00 GHz. The maximum 
time delay between 32.00 and 2.64 GHz is around 100 days for 
both flares, considerably longer than the estimated delays (be¬ 
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tween 20 and 70 days) for the other detected events. This might 
be an indication that the spectral evolution of the flares is tem¬ 
porarily modified by the onset of the main outburst. 

The largest variability brightness temperatures are measured 
at 2.64 GHz. They all exceed 2 x 10 12 K, implying a minimum 
Doppler boosting factor of ~ 4.3. In Fig. IB.ll we plot the time 
delays, the variability brightness temperature and the amplitude 
parameters for each flare separately. 

5.1.2. Flare decomposition for J0948+0022 

For J0948+0022, six prominent outbursts were identified in the 
light curves (Fig. [3] Table IB.21) . making the source a proto¬ 
type for applying the discussed method. As can be seen in 
Fig.EU for five of the six events the largest amplitude is seen 
at 32.00 GHz, while for flare 3 it is seen at 14.60 GHz. This is 
the result of high-frequency components that frequently appear 
at higher frequencies and expand fast towards lower bands. This 
behaviour is in accordance with the spectral evolution and the 
variability. The amplitude of the flares at different frequencies is 
shown in the upper panel of Fig. ED Interestingly, 

1 . flares 1 and 2 are remarkably similar to each other, but very 
different from flares 4, 5 and 6. The latter three are relatively 
isolated events, while the former two occur at the peak of the 
outburst phase. It is likely that this is exactly the reason for 
the difference in the amplitudes. 

2. flare 3 is different from all the others, with a large excess of 
flux at 14.60 and 10.45 GHz. The indications for an unusual 
spectral evolution of this event may be the interpretation of 
this phenomenology. 

On the other hand the time delays, between different fre¬ 
quencies shown in the middle panel of Fig. IB. ll reveal a mildly 
different behaviour only for flare 5. Unfortunately, a gap in the 
2.64 GHz data approximately where the peak is expected pre¬ 
vents us from excluding that the real time delays have to be 
shifted by about 30 days. Similarly, flare 3 at 2.64 GHz and flare 
5 at 32.00 and 23.05 GHz cannot be well constrained because of 
inadequate sampling. 

The largest 7' var is measured at 14.6 GHz. In the lower 
panel of Fig. ED we show 7' vilr versus v for the different flares. 
Most light curves give values of T vdI roughly between 2 and 
13 x 10 12 K. The differences from flare to flare are small (drops 
in the flux density are partially compensated for by decreasing 
timescales), so that they cannot provide independent estimates 
of the variability brightness temperature. The flare models at dif¬ 
ferent wavelengths, however, are independent of each other and 
hence the relatively small spread of 7’ var at different frequencies 
can be regarded as significant. The inferred values of 7’ var im¬ 
ply variability Doppler factors of the order of 5 to 7. The most 
prominent event occurred at around MJD 56444 (flare 3) and 
gave a 7’ var of more than 10 13 K implying a D of more than 8.7. 

5.1.3. Flare decomposition for J1505+0326 

Finally, for J1505+0326 we detected four main events during the 
monitored period as shown in Table IB.21 They are all relatively 
fast and of moderate amplitude, more similar to the character¬ 
istics of J0324+3410 than to those of J0948+0022. The highest 
variability brightness temperatures reach values of ~ 2 x 10 13 K 
and are detected at 2.64 GHz, implying Doppler factors of al¬ 
most 10.4. The spectral evolution of the flares seems to follow 
a standard pattern, with a long delay («;70 days) between 32.00 


Table 9. Variability parameters 7' var and D for the three sources 
whose datasets were long enough, and the frequency at which 
they are located. Here we report the highest values of 7' var . 


Source 

V 

(GHz) 

7 ’var 

(K) 

D 

J0324+3410 

2.64 

25 x 10 11 

4.3 

J0948+0022 

14.60 

13 x 10 12 

8.7 

J1505+0326 

2.64 

26 x 10 12 

10.4 


and 14.60 GHz and a considerably shorter one («40 days) be¬ 
tween 14.60 and 2.64 GHz. 


6. Radio spectral energy distributions 

Figure [4] shows the densely sampled broadband radio SEDs. 
There only data points with an S/N> 3 are shown. For guid¬ 
ance to the eye, the data points are connected with segments, and 
to avoid crowding the plots, only every fifth SED is colourecQ 
In Table [TO] a summary of their most characteristic parame¬ 
ters is given. These are the total number of available SEDs 
of at least two different frequencies, N; the average cadence 
quantified as the mean time-separation between two consecutive 
SEDs, mflA ( MJD A mm ( MJD > • qje coherence,that is, the mean time- 
separation between the earliest and latest measurement in a SED, 
max(MJD) - /win(MJD); and the spectral index in three sub¬ 
bands computed from the mean flux densities in those bands. 
The spectral index is defined as a — dlnS/d\n v, where S is the 
measured flux density at a frequency v. As can be seen there, the 
mean cadence (excluding J0849+5108) is shorter than a month, 
and the combined Effelsberg-IRAM spectra are made of mea¬ 
surements synchronous within less than two days. It is reason¬ 
able to assume that the shapes of the individual SEDs are free of 
source variability effects. 

From the plots in Fig. [4] it is evident that the mere sam¬ 
ple of these four sources displays th e main classes of spec - 
tral variability patterns described in lAngelakis et al.l (1201 2bl) . 
There, the variability patterns were interpreted as the result 
of the same underlying system consisting of a steep spectrum 
component attributed to a large-scale jet populated with high- 
frequency synchrotron comp onents that evolve acco rding to 
iMarscher & Gear! (Il985l) . For lAngelakis et akl (1201 2 hi) the dif¬ 
ferent observed classes were governed by the amount of “activ¬ 
ity” that the observing bandpass allows us to see, and was also 
caused by different intrinsic properties. Except for J1505+0326, 
where the spectral evolution is not discernible, all other three 
NLSyl would be classified as type 1 or 2 in that classification 
scheme because they show an intense spectral evolution with or 
without a steep - presumably quiescent state - spectrum. 

6.1. J0324+3410 

For J0324+3410 Fig.|4]clearly shows the emergence of two dis¬ 
tinct components: 

- A steep spectrum component that appears at frequencies be¬ 
low 10.45 GHz and mostly retains a constant negative spec¬ 
tral index that reaches values as low as almost ~ -0.5 (see 

1 Animated SEDs can be found at 

http://wwwB.mpifr-bonn.mpg.de/div/vlbi/fgamma/ 
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Table 10. Mean characteristics of the observed radio SEDs shown in Fig. [4] The spectral indices are computed in three sub-bands, 
and along with the mean values, the extremes are reported as well. 


Source N Cadence Coherence low: 2.64 - 8.35 GHz mid: 10.45 - 23 GHz high: >43.05 GHz 




(days) 

(days) 

min 

mean 

max 

min 

mean 

max 

min 

mean 

max 

J0324+3410 

55 

25 

2.0 

-0.49 

-0.21 

+0.12 

-0.53 

-0.04 

+ 1.57 

-0.67 

+0.07 

+0.77 

J0849+5108 

12 

84 

0.02 

+0.04 

+0.18 

+0.44 

+0.08 

+0.32 

+0.66 

+0.28 

+0.63 

+ 1.08 

J0948+0022 

75 

26 

2.2 

-0.13 

+0.51 

+ 1.26 

-0.83 

+0.27 

+ 1.02 

-1.73 

-0.03 

+ 1.58 

J1505+0326 

43 

29 

0.2 

-0.22 

+0.13 

+0.34 

-0.66 

-0.18 

+0.18 

-0.86 

-0.24 

+0.07 


Notes. Column description: 2: the total number of available SEDs; 3: the mean cadence i.e. the separation between two consecutive spectra 
measured in days; 4: is the max(MJD) - mm(MJD) within one single SEDs; 5, 6, 7: are the minimum, the mean and the maximum spectral indices 
for the low band and so forth. 


Table [JOj and could be the optically thin part of a syn¬ 
chrotron self-absorbed (SSA) component. 

- A sequence of high-frequency components (HFC) that intro¬ 
duce intense variability as a result of their spectral evolution. 

A convenient way to understand the reason for the observed 
behaviour is to examine the activity in terms of spectral index 
distribution in three different frequency sub-bands: 

Fow: for v e [ 2.64 , 8.35] GHz 
Middle: for v e [10.45 , 23.05] GHz 
High: for v e [32.00,142.33] GHz. 

In Fig. [5] we show the spectral index distributions in these 
sub-bands. The histograms are normalised so that they comprise 
the probability density functions (PDF) of the three indices. 
The bin size was selected to be similar (if not equal) to the 
mean standard deviation computed over the corresponding 
dataset. This representation compresses the behaviour and the 
phenomenology discussed earlier and immediately indicates the 
sub-band where the activity is mostly concentrated. 

As can be seen in Fig. [4] the region above roughly 10.45 GHz 
hosts all the source activity in terms of spectral evolution. It ap¬ 
pears to be intensely variable, with the outbursting events emerg¬ 
ing at the upper end of the band and transversing downwards 
following a characteristic evolutionary path; first a build-up at 
high frequencies, progressively harden as the flux increases, and 
a descent towards the left end of the bandpass while shifting the 
peak to lower flux densities to eventually diffuse in the persistent 
steep spectrum below 10.45 GHz. 

This scheme agrees at least qualitatively with the evolution¬ 
ary scenario proposed bv lMarscher & Gelul G 985 ). The succes¬ 
sion of these components and the superposition of their evolving 
spectra results in a flat average spectrum - with a spectral index 
practically around 0.0 - which at times can be as hard as +1.6 
(see Table ITOli. 

The steep spectrum component that dominates the SEDs at 
lower frequencies is not static. It shows variability in a self¬ 
similar fashion and mostly retains a steep character. The spectral 
index is positive for less than 7% of the time. The lowest in¬ 
dices (« -0.5, see also Fig.0 indi cate the canonical values seen 
in large-scale relaxed jets (e.g. [Angelakis et all 120091) . A natu¬ 
ral way to explain this behaviour is to assume that below about 
10 GHz, the high-frequency components are optically thin. They 
transverse that part of the bandpass and slowly fade away and 
gradually contribute less and less power. This lets the overall 
spectrum appear to move downwards in flux density while pre¬ 
serving a rather constant index. This interpretation implies that 
high-frequency components with similar intrinsic parameters are 


generated (i.e. magnetic field and particle density). At the very 
low frequencies (2.64 and 4.85 GHz), the moderate spectral vari¬ 
ability indicates an underlying large-scale jet. Essential in favour 
of this scenario would be the detection of significant polarisation 
in that part of the emission, as we discuss in Sect. [8] 

The upper end sub-band (v > 43.05 GHz) can become as 
steep as of a ~ -0.7 (Fig. [5}. denoting that the transition from 
optically thick to thin regime of the HFC occurs in that spectral 
region. This occurs rather rarely; only 14% of the cases show 
such a steep index, most likely because as an event reveals its 
optically thin part, a new one emerges, that fast hardens the oth¬ 
erwise softening spectrum. 

In the middle sub-band (10.45 < v < 23 GHz), however 
the source spends half the time in soft and half in hard state 
since the “mixing” of thin and thick events is prominent there. 
In any case, the transition from thick to thin has very serious 
consequences because if proven true, it should be accompanied 
with exactly 90°rotations of the projected electric vector position 
angle (EVPA) as is described bv lMvserlis et al.l (120144) . 

6.2. J0849+5108 

The limited dataset available for J0849+5108 (see Fig. [4} al¬ 
ready reveals a very interesting behaviour. First to be noted is 
the complete absence of any steep spectrum signature from the 
SEDs that have been gathered so far. The spectrum oscillates 
between flat and hard within all three sub-bands. As an exam¬ 
ple, the low sub-band spectral index varies between +0.02 and 
+0.44. In the high sub-band it reaches values higher than +1. 
Furthermore, the events always show their optically thick part, 
which - in the framework of incoherent synchrotron processes - 
is expected to be associated with an insignificant degree of lin¬ 
ear polarisation. Clearly, more duty cycles are needed to allow 
better quantification of its spectral behaviour. 

6.3. J0948+0022 

The complete absence of a quiescence spectrum also holds for 
J0948+002. The spectra appear to be highly inverted and in¬ 
tensely variable. An impressive spectral evolution with HFCs 
emerging in the high sub-band is observed. The HFCs often 
transit to optically thin phases within the middle sub-band. 
Interestingly, the lowest band spectral index (2.64 < v < 
8.35 GHz) oscillates between a minimum of -0.13 and a max¬ 
imum of +1.26, while the middle- and high-band indices oscil¬ 
late between -0.83 and +1.02 and -1.73 and +1.58, respectively 
(c.fFig ,[6}. This implies that we witnes the spectral evolution of 
the HFCs, but fail to detect the quiescent steep spectrum, if any. 
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Fig. 5. Probability density functions of the spectral indices ob¬ 
served for J0324+3410 in three different sub-bands: low, mid¬ 
dle, and high. The bin size is selected to be approximately one 
standard deviation of the plotted indices. 
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Fig. 6. Probability density functions of the spectral indices ob¬ 
served for J0948+0022 in three different sub-bands: low, mid¬ 
dle, and high. The bin size is selected to be approximately one 
standard deviation of the plotted indices. 
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Fig. 4. Radio SEDs of the four monitored NLSyl galaxies. The 
data points are connected with straight segments to guide the 
eye. The legend denotes the fractional year for each SED. 


VLBI studies (e.g. Karamanavis et al. in prep.) also reveal only 
a compact core structure in agreement with the previous finding. 

Furthermore, the relatively fast pace at which the spectral 
evolution occurs is interesting and rather exceptional. In fact, 
significant variability has previously been detected with a mean 
cadence of 25 days. This is also seen in the light curves, as dis¬ 
cussed in Sect. 14.31 Figure [6] again gives the probability density 
functions of the three spectral indices as for J0324+3410. 

6.4. J1505+0326 

The behaviour of J1505+0326 is somewhat different from the 
behaviour of the sources discussed previously. Most of the time, 
it shows a convex spectrum that has a global maximum in the 
vicinity of 5 - 10 GHz and occasionally seems to be modulated 
by the evolution of HFCs that become optically thin. For that 
rea son, it could be classified as type 5b in the scheme proposed 
by lAngelakis et al. ( 2012c ). The low sub-band spectral index 
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Table 11. Radiative and kinetic jet powers calculated from the 
14.60 GHz observations averaged over the period discussed 
here. Literature values of their black-hole masses are as well. 


Source 

Rj,rad 

(erg s' 1 ) 

P j,kin 

(erg s' 1 ) 

P j,tot 

(erg s' 1 ) 

logM 1 
(log M 0 ) 

J0324+3410 

1.7 x 10 43 

3.0 x 10 43 

4.7 x 10 43 

7.6 

J0849+5108 

6.6 x 10 44 

2.5 x 10 45 

3.1 x 10 45 

7.5 

J0948+0022 

7.9 x 10 44 

3.0 x 10 45 

3.7 x 10 45 

7.9 

J1505+0326 

5.0 x 10 44 

1.7 x 10 45 

2.2 x 10 45 

7.3 


References. (l)[Foschini et~aTI 120141 ) 


reaches moderately hard states (max (a) = +0.34), while the 
high sub-band can become very soft (min (a) = -0.86). 


7. Jet powers 

Here we wish to quantify the power delivered by the four 
NLSyls in our sample in the radio bands and examine whether 
this could be delivered by a jet similar to those in typical blazars. 

Within the jet theor etical framework proposed by 
lBlandford & Konigl (1979 1 and empirical relations that 
have been proposed for relativistic iets in AGNs and galactic 
binaries bv IFoschini (i 201 la . 2012a . 2014), the radio measure- 
ments_can accomm odate the computation of the jet power (see 
also iFoschini et al.ll2014l) . In fact, it is possible to calculate both 
the radiative and kinetic (electron, proton, magnetic) jet power 
from the radio core emission. In a conical jet, the optically thick 
radio emission is linked to the jet power as 

Tradio oc Fj 7/12 , (3) 

By studying a sample _of 88 y-ray AGN (FSRQs, BL Lac 
Objects, RLNLSlsl. lFoschi~ni](20l4 ) found that the best normal¬ 
isation coefficients values are 

log Fj.rad = (12 + 2) + (0.75 ± 0.04) log L rad i 0 ,core (4) 

log 

^j,kin — (6 ±2) + (0.90 ±0.04) log U 

adio,core ? (5) 

where P h ai and Pj.km are the radiative and kinetic (proton, elec¬ 
tron, and magnetic field) powers of the jet, while L ra dio,core is 
the A'-corrected luminosity of the radio core at 14.60 GHz. The 
A'-correction was made by multiplying the monochromatic flux 
by the factor of (1 + z) a ~ 1 , assuming a flat radio spectral index 
(a ~ 0, see also Sect. 01. 

Although it would be possible to use different frequencies 
as well, we persist in using 14.60 GHz for a number of rea¬ 
sons, such as to be consistent with the measurements used to 
extract the used normalisation coefficients; the good sampling of 
the corresponding light curve; the fact that this band is less con¬ 
taminated by extended emission, which gives more reliable esti¬ 
mates of the radio core emi ssion; it allows comparison with other 
proje cts such as MOJAVE (iKellermann et all 120041 : [Lister et all 
l2009ll2()Tl . 

In Table |TT] we provide the naturally weighted averages of 
the radiative / J j ra d, kinetic / J ]km , and total powers Pj, to t- The av¬ 
eraging was made in the power rather than the flux density do¬ 
main. In the same table we present the most recent estimates for 
the black-hole masses (Foschini et all [2014 1. 


7.1. J0324+3410 

The prominent event seen in the source light-curves at around 
MJD 56300 (January 201 3) gave a co rr espon ding increase in the 
radio power at that time. iTibolla et al.l (2013b reported a change 
in the X-ray spectrum from January to February 2013 associated 
with that event, with the emergence of a hard tail during the sec¬ 
ond epoch. The average values of the radiative Pj, ra d and kinetic 
jet power / J | km - as measured in the present work - are 1.7 x 10 43 
and 3.0 x 10 43 erg s 1 . Th e sam e quantities calculated from the 
SED modelling bv lAbdo et al.l (2009d) after the first year of the 
satellite operations were 6.3 x 10 42 and 2.2 x 10 44 erg s' 1 . 


7.2. J0849+5108 


Unfortunately, the limited dataset available for J0849+5108 
is not adequate to reveal the source activity characteristics. 
Furthermore, available measurements are no t simul taneo us t o 
the y-ray active phase it went through in 2010 ( Foschinil 201 lbh. 
The radiative jet power reported by iD’Ammandoet al.l (20121 
l2013cb is around (5 - 50) x 10 44 erg s 1 (the maximum refers to 
the outbursts), which is consistent with our calculations, which 
gave an average of 6.6 X 10 44 erg s 1 . 


7.3. J0948+0022 

J0948+0022 was extensively studied bv iFoschini et al.l (2012b . 
who discussed the results of more than three years of multi¬ 
wavelength data over the period between 2008 and 2011. They 
presented the calculations of jet powers in terms of SED mod¬ 
elling. With the current work - which comprises the natural ex¬ 
tension of that monitoring to mid-2014 - we can compare these 
findings with those of the currently applied method for the over¬ 
lapping period. A caveat of the SED modelling has been the 
delay of the radio emission with r espect to the optical-to-y-ray 
emission, as is generally expected. iFuhrmann et al . (12014 ) found 
that the average time delays at the source rest frame can span 
from less than 10 to more than 70 days, systematically decreas¬ 
ing towards higher frequencies, which provides yet another ar¬ 
gument in support of using the 14.60 GHz for the current calcu¬ 
lations. 

The radiative powers computed by IFoschini et al.l (2012b 
were based on SEDs that were selected on the basis of the y-ray 
activity. During the lowest activity SED, the jet radiative power 
was estimated to be Pj, ra d = 0.4 x 10 45 erg s' 1 , while during the 
highest activity state it was Pj, ra d = 7.9 x 10 45 erg s _1 . From 
our radio dataset, the radiative power ranges between 0.4 and 
1.8 x 10 45 er g s 1 . For the jet kineti c pow er, the SED modelling 
applied by iGhisellini & Tavecchiol (2009b showed that it ranges 
from 0.2 to 4.1 x 10 47 erg s' 1 , while based on the 14.60 GHz 
measurements of the present work, it varies in the range be¬ 
tween 1.5 to 7.8 x 10 45 erg s' 1 . This result is very interesting 
when the fundamental assumptions are considered. The kinetic 
jet power - as calculated from the SED modelling - relies on the 
assu mption that there is one p roton associated with each elec¬ 
tron ( Sikora & Madeiskl . 12000 ). If that were not the case and in¬ 
stead there were one proton for more electrons, the kinetic power 
would be much lower, which is what our finding indicates. 


7.4. J1505+0326 

The multi-frequency light curves of J1505+0326 show at least 
one major event with mild spectral evolution. Despite its mild 
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variability, we can c ompare, o ur values with those from SED 
modelling reported in lAbdo et al.l d2009ch . 1.0 x 10 44 and 1.6 x 
10 46 erg s -1 are the values of the radiative and kinetic jet powers. 
The averages we calculated are 5.0 x 10 44 and 1.7 X 10 45 erg s _1 , 
respectively. The almost constant activity of the source is con - 
sistent with what was reported by ID’Ammando et all d2013ai) . 
who analysed the monitoring at 15 GHz performed by the Owens 
Valley Radio Observatory (OVRO). 


8. Polarisation 

Assuming that the dominant radiation mechanism is incoherent 
synchrotron processes and that the variability is attributed to 
evolving synchrotron self-absorbed components, it is expected 
that significant linear polarisation is detected. In this context, the 
degree of polarisation and the EVPA would depend on the emit¬ 
ting pl asma opacity and hence the spectral index dPacholczvkl 
il970L 119801). Assuming high ly ordered and intense magnetic 
fields dSciama & Reeslll967h . as those expected during flares, 
even circular polarisation could occur. Furthermore, the orienta¬ 
tion of the electric vector should also depend both on the opti¬ 
cal thickness of the material and on the assumed magnetic field 
configuration. Multi-band polarisation studies are then a unique 
probe of the underlying physical conditions at the emission re¬ 
gion. 

For the sample discussed here, iNeumann et al.l (11994 1 
found optical p o larisat ion from J0324+341 0, while 
iMoore & Stockman (i!981h and ISitko et all (1 19841) detected 
“high and variable optical polarisation” from J0849+5108. This 
has been seen as an indi catio n th at a relativistic jet is operating at 
the source. Similarly. iDoi et ak| (12012! ) have detected significant 
polarisation in the radio. 


8.1. Radio polarisation 

lAngelakis et al.l (i2013l) explored the possibility of detecting 
radio linear polarisation at the intermediate frequency of 
8.35 GHz. The data in Table 3 there were corrected for in¬ 
strumental polarisation and indicated that all four sources dis¬ 
play detectable polarisation of a remarkably low magnitude. 
Exceptionally, J0324+3410 showed values higher than those 
typically observed in A GN a t these frequency bands (~ 4 % at 
5 GHz, e.g. lKlein et al.l2003l) . Its fractional linear polarisation at 
this frequency is about 6 % with an EVPA of about 36°, while at 
4.85 GHz it reaches values of 7 % with an EVPA of about 44°. 
This phenomenology would agree with assuming that the opti¬ 
cally thin emission dominates towards lower frequencies. The 
4.85 GHz emission is less contaminated by the optically thick 
part of an HFC and displays a higher fractional polarisation. At 
both frequencies the EVPA appears to be almost perpendicular 
to the jet orientation (the jet axis is about 120°, Karamanavis et 
al. in prep.), so that the projected magnetic field appears to lie 
almost parallel to the jet, possibly implying a helical jet with a 
very long helix step. 

Similar arguments hold for J0849+5108, although the effect 
is much more moderate. Unlike J0324+3410, the other NLSyls 
show a very low polarisation degree. As is indicated by their 
spectral variability pattern, the sources undergo flaring events 
that show their optically thick part at 8.35 GHz. 


8.2. R-band polarisation 

First attempts to monitor the R-band linear polarisation in 2013 
with the RoboPol instrument (iKing et all 120141 iPavlidou et all 
120 1 4l) showed that J0324+3410 and J0948+0022 present a rather 
low fractional polar isati on that does not exceed a few percent 
(=s 3 %). Iltoh et all (.120141 ) found similar values for the former 
case, while the latter showed past events_ of polarisation values 
that reached almost 40 % (Iltoh et al.Ll20H . For J0849+5108 we 
detected polarisation levels of about 10 %. 

9. Discussion 

The main goal for the analysis has been to study the relativis¬ 
tic jets that emit the observed radio power and to quantify its 
properties. The discussed dataset comprises the longest term 
multi-wavelength radio monitoring datasets of the known y-ray- 
detected NLSyls to date. 

The focus was placed mostly on the variability proper¬ 
ties of the radio emission. In this context, we follow different 
paths. We introduced a flare decomposition algorithm to quan¬ 
tify the variability amplitude, the involved timescales, and its 
frequency dependence for each flare separately. We then ex¬ 
amined whether variability retains the same characteristics over 
source, frequency or event. Subsequently, the variability param¬ 
eters (amplitude and time scale) served as the base for comput¬ 
ing the limiting values for the variability brightness tempera¬ 
tures and the corresponding Doppler factors. In the frequency 
domain, much attention was given to the spectral evolution seen 
in each source, chiefly because of the signatures that the vari¬ 
ability mechanisms leave in the evolutionary path followed by 
the flaring events. Furthermore, the power delivered by source 
in the radio regime was computed at one selected frequency and 
was compared with the powers seen in the jets of typical blazars. 
Finally, the radio and the optical polarisation was investigated as 
a signature of relativistic jet emission. 

The mean flux densities are between 200 and 600 mJy. All 
sources appear significantly variable at cm and mm bands, ex¬ 
cept maybe J0849+5108, for which it is rather premature for 
such a statement to be made. The flare decomposition algorithm 
revealed several events in each light curve. The characteristics 
of different events can vary even in the same light curve, im¬ 
plying that either the underlying mechanism is not the same, or 
that each time different evolutionary stages are detected. The de¬ 
lays between the peaks at different frequencies of the same event 
range between a few days, tenths of days, and 180 days, showing 
a power-law frequency dependence. The flare timescales range 
between ~ 30 and 200-300 days. The computed highest vari¬ 
ability brightness temperature 7' var is between 10 12 and 10 13 K. 
Consequently, values ranging from 4.3 to 10.4 were found for 
the corresponding Doppler factor, D. 

From examining the behaviour in the frequency domain, it 
is immediately clear that intense variability occurs in the same 
ways as routinely seen in blazars. Specifically, the variability 
events seem to appear at high frequencies and evolve downwards 
in frequencies, following an evolutionary path seen in blazar 
flares. Interestingly, the power of the events is a significant part 
of the overall power. Apart form the case of J0324+3410 that 
shows signs of a steep underlying spectrum the rest give no hints 
for such in our bandpass. 

The 14.60 GHz radio power ranges roughly between 10 43 
and 10 45 erg s in agreement with what should be expected 
from a relativistic jet. The 8.36 GHz polarisation data indicated 
low degree of fractional polarisation for three of our sources. For 
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J0324+3410 however it reached values up to almost 6 %. Strong 
indication that we are dealing with the non-thermal emission 
originating at a relativistic jet which is even dominating over 
a core emission despite the source compactness at VLBI scales. 
Despite the interesting findings of first attempts to conduct op¬ 
tical linear and radio linear and circular polarisation study (that 
we have already discussed in Sect. [8] and will discuss below), a 
methodical and complete polarisation analysis will be presented 
in subsequent publication after longer optical polarisation time 
baselines are completed. 


9.2. J0849+5108 

Unfortunately, the limited number of observations that have been 
performed on the source do not allow a proper light curve anal¬ 
ysis. The radio SED, however, already reveals clear signs of - 
at least mild - spectral evolution indicative of recursive activity, 
a typical characteristic of radio blazars. Its spectral index shows 
variability at all sub-bands (see Sect. [6]) and reaches values of 
+ 1, indicating that our bandpass samples a part of the spectrum 
where most of the radio activity occurs. The dominance of the 
assumed optically thick emission could also explain the very low 
radio polarisation. 


9.1. J0324+ 3410 


As we extensively discussed in Sect. [4] the source light-curve 
has similar characteristics as blazars with the only exception of 
the moderate average flux densities. In combination with the pic¬ 
ture drawn from the radio - almost monthly sampled - SEDs 
and their measured spectral indices, this points to the scenario 
of a jet within which high-frequency components are generated. 
The HFCs subsequently evolve transversing the possible opac¬ 
ity phases, b efore disappearing with an expansion toward s low 
frequencies dMarscher & Gear! 119851 : iTiirler et~ail. 120001) . The 
emergence of later events of similar phenomenology indicates 
that the physical conditions are stable (particle densities, mag¬ 
netic field strength) at the location of the emission and could be 
in accordance with a reoccurring activity that would result from 
a stable or, a slowly evolving acceleration region. The variability 
of J0324+3410 is characterised by fast variations - shorter than 
65 days, even at the lowest frequencies - with relatively low 
amplitude. From the flare decomposition exercise we estimate 
that the highest brightness temperature limits are all higher than 
2 x 10 12 K, implying a Doppler boosting factor of at least 4.3, 
rather lo w com pared with the variability values seen in blazars 
(e.g. lHovatta et al.Ll2009h . 


An essential element of the exact processes operating at the 
source would be studying of polarisation behaviour at radio and 
optical wavelengths. As was discussed earlier, J0324+3410 dis¬ 
played fractional linear polarisation o f about 6 % and 7 % a t 
8.35 GHz and 4.85 GHz, respectively dAngelakis et al.Ll2012a ). 
This increase towards lower frequencies might be an indica¬ 
tion that optically thin emission - arguably from a relativistic 
jet - becomes progressively dominant towards lower frequen¬ 
cies. At higher frequencies the emission is dominated by op¬ 
tically thick (and hence less polarised) emission of an evolv¬ 
ing HFC. In the optical regime the detected significant polar¬ 
isation has been seen as ev idence of an, at least mildly, rela- 
tivistic jet ( Yuan et al.ll200 8l). Our R-ban d RoboPol monitoring 
dKing et al. . 120141 : Pavlidou et all 120141) showed that the opti¬ 
cal polarisation of the source is very low and shows fluctuations 
between 1 and 3 %. Its /Gbaud magnitude is almost stable at 
roughly 15.3 mag. A possible explanation for this unusually low 
value of the optical polarisation may be the contamination of 
the emission with an unpolarised stellar component that would 
result a lower net polarisation. 


Finally, the total jet powers computed from our 14.6 GHz 
dataset agree well with those computed with SED modelling that 
accounts for all different emission components. In conclusion, 
there is compelling evidence that the source radio emission is 
dominated by a jet emission that shows characteristics often seen 
in blazars, but not necessarily as powerful. 


9.3. J0948+0022 

The detection of J0948+0022 as a y-ray emitter (lAbdo et all 
l2009al) triggered the monitoring programme we discussed here. 
Consequently, its light curves are the best sampled. It is very 
radio-loud (iZhou et all 12003 ). in contrast to the typical NLSyls, 
and displays light curves with clear intense variability. Its radio 
SED shown in Fig. [4] reveals a very interesting case with se¬ 
quential HFC appearing at high frequencies, evolving through 
the bandpass before fading with the emergence of a new com¬ 
ponent. The lowest band spectral index shows values as inverted 
as +1.26, while the high band spectral index can be as steep as 
-1.73, implying that the source activity occurs within our band¬ 
pass. A careful cross-correlation of spectral index trains with the 
associated radio polarisation will be key to really proving the 
presence of the HFC that cause the variability; this would favour 
a typical blazar-like jet. 

From the flare decomposition it appears that the best con¬ 
straint in the variability brightness temperature is set by the 
14.60 GHz measurements and gives values 13 x 10 12 , imply¬ 
ing a variability Doppler factor of ~ 8.7. Karamanavis et al. 
(priv. comm.) have been trying to detect moving components 
in 15 GHz MOJAVE data for the purpose of measuring appar¬ 
ent speeds. However, the source appears to be unresolved. In 
conclusion, it could be said that this source resembles most the 
typical blazar behaviour of the four sources in our sample. 

In Sect. I7.3l we discussed that the jet powers computed from 
SED modelling dFoschini et all 12012i) are more than an order 
of magnitude larger that those computed here. Interestingly, this 
could be an indication that the jet composition is different from 
the compositi on assumed for SED mode lling. 

Recently, ID’Ammando et al] (1201 4l ) reported the detection 
by Fermi of a strong y-ray flare observed between December 
2012 and January 2013 in the 0.1-100 GeV energy range. Its 
peak occurred at MJD 56293. A quasi-simultaneous flare was 
observed from optical to X-rays, suggesting a tight connection of 
the corresponding emitting regions. The examination of OVRO 
15 GHz light curves showed no flaring activity at the time of the 
y-ray flare. A strong radio peak was instead observed in June 
2013. The authors hypothesised that there might be a delay be¬ 
tween the y-ray and the radio activity of about five months, due 
to opacity effects and propagation of the shock along the jet. The 
case becomes very interesting when examining our light curves, 
which show the detection of a radio-flare quasi-simultaneous to 
the y-ray flare. As can be seen in Fig. [4] the 23.05 GHz data show 
a very prominent flare. It starts at the beginning of December 
2012, peaks between December 15 and January 4, 2013, and 
ends at the beginning of February 2013. The flare is weaker in 
the 32.00 GHz light curve, while at 14.60 GHz it is practically 
invisible. Unfortunately, the insufficient sampling of highest fre¬ 
quencies in our programme prevent us from quantifying the flare 
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dependence in frequency. Our data suggest that the y —ray flare 
propagates to radio frequencies down to ~23.05 GHz with al¬ 
most no time delay, showing that the propagation of the shock 
along the jet is fast. At lower frequencies, the flare drops dras¬ 
tically, probably because of opacity effects. The examination of 
the radio SED implies a subtle HFC that barely reaches the cen¬ 
tre of our bandpass before it disappears. In any case, the sudden 
disappearance of the event as a function of frequency demands 
more investigations of the mechanisms behind the very broad¬ 
band emission (from y rays to radio) and its evolution. 

9.4. J1505+0326 

The radio SED of J1505+0326 differs from the usual behaviour 
seen in the other NLSyls in our sample. Its spectrum is most of¬ 
ten convex with its peak in the range of 5-10 GHz. Occasionally, 
there appear signs of evolution caused by an HFC entering the 
bandpass. It is essential, however, that the total number of avail¬ 
able cycles does not allow deciding whether the spectrum is in¬ 
trinsically convex or whether it possibly changes at a slow pace. 

The flare decomposition detected four rather short events and 
have moderate amplitude despite the source’s radio-loudness. 
The highest variability brightness temperature limit was com¬ 
puted at 2.64 GHz and yielded values of ~ 2.6 x 10 13 K, imply¬ 
ing a boosted relativistic jet with Doppler factors of more than 
10 . 

10. Conclusions 

Before we conclude, we emphasise that our findings cannot nec¬ 
essarily be assumed to be representative of all RL, y-ray NLSyl 
galaxies. The most outstanding findings of this work are as fol¬ 
lows: 

1. We presented the most systematic radio monitoring of RL 
and y-ray-detected NLSyls to date. The obtained light 
curves spread over periods longer than five years. The ac¬ 
quired dataset includes eight frequency bands between 2.64 
and 142.33 GHZ and is available online for the community 
for further studies. 

2. The light curves show the typical phenomenology seen in 
blazars, although they generally have lower flux densities. 
All of them show variability events that in most cases have a 
smaller amplitude, however. Sometimes very prominent and 
energetic events appear and dominate the light-curve appear¬ 
ance. 

3. We introduced an alternative method for detecting flares and 
quantifying their characteristics, which - contrary to the tra¬ 
ditional methods - is approaching the problem on the basis of 
individual flares. In that context, the properties of each flare 
can change over frequency; similarly, the different events in 
the same light curve can appear with different properties, 
which indicates either different variability mechanisms or a 
different flare evolution stage. 

4. The flare decomposition method showed that the computed 
variability brightness temperatures are in general moderate 
- in contrast to the majority of blazars -, implying Doppler 
factors that do not exceed about 10. The operating jet then 
must be moderately relativistic. 

5. The same method showed that the duration of the events is 
somewhat shorter than what is seen in blazars, resulting in 
slightly more frequent events. 

6 . For the three cases where longer periods have been covered, 
the radio SEDs clearly show intense variability that is mostly 


characterised by clear spectral evolution of the events. The 
evolution seems to occur slightly faster than in blazars - as 
was also indicated by the results of the flare decomposition 
- but follows patterns possibly caused by evolving internal 
shocks. Furthermore, in one case the shape of the radio SED 
towards the lower frequencies gives a direct indication of a 
quiescent jet (steep spectrum, steady). 

7. The jet power estimations gave values in accordance with 
jet outputs comparable to the least energetic blazars (BLLac 
objects). For 10948+0022 the comparison with SED mod¬ 
elling results indicates lower values that may be attributed to 
different p-e ratios in the jet. 

8 . The first multi-frequency radio polarisation results showed 
that our sources have shown very low or negligible polarisa¬ 
tion. In contrast, J0324+3410 a displayed polarisation degree 
of 6 and 7 % at 8.35 and 4.85 GHz. This is rather high, es¬ 
pecially given the compactness of the source, and could be 
explained by the jet dominating the core in the total power. 
The increase towards lower bands is in accordance with that 
assumption as more and more optically thin emission from 
an underlying jet becomes detectable. From the above it be¬ 
comes rather evident that 

9. a mildly relativistic boosted jet must be operating at these 
sources and probably is responsible for the observed radio 
emission. This jet lie, probably towards the low end of the 
energy distribution, however. 

10. The characteristics of its radio emission are qualitatively 
similar to those in blazars, but are generally less intense. The 
spectral evolution seems to occur at a faster pace, involv¬ 
ing shorter timescales probably related to the systematically 
smaller black-hole masses or higher accretion rates. 

11. In the particular case of J0948+0022, a flare detected by 
Fermi towards the end of 2012 seems to be present in our 
23.05 GHz dataset, but is practically absent at frequencies 
below that. This unusual behaviour deserves further investi¬ 
gation of the possible mechanism that could cause such an 
abrupt flare disappearance. 
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Appendix A: flare decomposition method 

The detection and parametrisation of different flares (which are 
then used to compute the associated variability brightness tem¬ 
perature) requires solving several problems, most of which have 
to do with the paradox that before an independent parameteri- 
sation of all flares is reached, we need to identify their average 
shape and spectral evolution, minimise the arbitrariness of the 
results. The problems to be addressed can be summarised as fol¬ 
lows: 

1. We ned to identify of the strongest flare, accounting for the 
fact that its amplitude may strongly depend on frequency so 
that the strongest flare at one frequency may be a minor event 
at others. 

2. We also need to localise of the same flare at different fre¬ 
quencies, given its spectral evolution. 

3. The flares need to be modelled efficiently using the minimum 
number of free parameters. 

4. Finally, we have to test and possibly correct the model to 
optimise the results. 

The implementation of these steps was made as follows: 

1. We estimated the average time delays between the flux den¬ 
sity variations detected in different light curves. The main 
complication in this procedure is t he fact that using a stan- 
dard cross-correlation function (see lEdelson & Kroliklll988t 
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iLehar et all il992 i on pairs of light curves does not ensure 
consistent results. For example, assuming the estimated de¬ 
lays between three light curves a, b, and c to be r a ,b, T u, c 
and r b, c , it might well be that r aiC 4- T a ,b + t^ c - A solution 
of this problem is to include all the possible pairs of light 
curves when calculating the time delays. By visual inspec¬ 
tion, a maximum time delay between light curves, T max , is 
identified. Then, for every pair of light curves i and j with 
i < j, a cross-correlation function r,j(T;j) is calculated for 
any time delay ry in the range (-r max , T max ). For n moni¬ 
tored frequencies an overall correlation function can be cal¬ 
culated as 

t 2 ,3,.... = Y\ r u( T /j)» (A. 1) 

i,j 

where 

Ti,j = Ti, 1+1 + Ti+ij+2 + ... + Tj-1J, 

We computed R for any possible set of time delays, exclud¬ 
ing all the combinations involving < 0. The most 

significant set of time delays was selected as the one provid¬ 
ing the highest overall correlation value. 

2. We composed a “cumulative” light curve, in which the 
strongest flare was identified. The light curves were shifted 
in time by the best fit delays found in the previous step, 
so that their variations were aligned. Their fluxes were nor¬ 
malised by their standard deviation, to have similar varia¬ 
tions at all wavelengths. Assuming for example that all the 
variations at 86.24 GHz have an amplitude of 1 Jy and at 
2.64 GHz have an amplitude of 0. Jy, they all count the same 
for the identification of the flares. This way, we ensured that 
our detection algorithm was not biased towards the high¬ 
est frequencies. The resulting light curves were merged to 
create one well-sampled total “cumulative” light curve. The 
merging of the light curves does not require averaging; the 
points from all the light curves are combined regardless of 
the frequency at which they were observed. The previous 
steps are meant to make all points equivalent, giving us the 
opportunity create a well-sampled light curve in which every 
data point from every frequency participates with the same 
weight. To construct the cumulative light curve, one could 
identify the following steps: 

3. We identified the strongest flare. The cumulative light curve 
was smoothed (with an adaptive smoothing factor) and lo¬ 
calised the highest flux density. 

4. We localised the same flare at different frequencies. 
Knowing the average time delays between variations occur¬ 
ring at different frequencies, the strongest flare can be iden¬ 
tified in each of the actual light curves. Gaussian fit mod¬ 
els, different for each wavelength, were applied to the actual 
light curves to determine the amplitude and timescale of the 
strongest flare; they were then subtracted from the original 
flux densities, creating a new set of light curves (one per 
wavelength) for which all the steps from step 3 on are re¬ 
peated. The procedure was iterated until no other significant 
flares are detected. 

5. Model test. The residual light curves were checked to verify 
that no peaks or dips are systematically present in the posi¬ 
tion of the subtracted flares. Their existence would indicate 
that the flares are not properly modelled; new models would 
be created and the whole procedure repeated until a satisfac¬ 
tory result is reached. 

This approach has the advantage of individually localising 

and modelling each event. Once a flare is detected in a cumula¬ 


tive light curve, the average time delays returned by step 1 al¬ 
low us to infer its approximate position in the real light curves. 
The accurate position is then identified from the flux peak clos¬ 
est to the approximate position. In this way, we are able to es¬ 
timate the actual time delays as a function of wavelength for 
each flare separately, and therefore detect possible variations in 
their spectral evolution. Naturally, the procedure also provides 
information about how the amplitude and the rising and decay¬ 
ing times of the flares vary with wavelength, which leads to a 
full parametrisation of the flux density variations observed in the 
sources; this is necessary to compute brightness temperatures. 

Appendix B: Results of the flare decomposition 

Here we append the results of the flare decomposition method 
described in Sect. EU 
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Table B.2. Results of the flare decomposition for all three sources. 


V 

(GHz) 

Time 

delay 

(d) 

Amplitude Peak 

position 
(Jy) (MJD) 

Time 

scale 

(d) 

Tv ar 

(10“ K) 

D 

32.00 

-100 

1.284 

J0324+3410 

Flare 1 
56313 

46 

4.6 

2.2 

23.05 

-101 

0.744 

56312 

38 

8.0 

2.7 

14.60 

-105 

0.514 

56308 

33 

18 

3.6 

10.45 

-82 

0.205 

56331 

47 

6.9 

2.6 

8.35 

-81 

0.133 

56333 

65 

3.9 

2.2 

4.85 

-52 

0.105 

56362 

64 

8.6 

2.9 

2.64 

0 

0.081 

56413 

61 

25 

4.3 

32.00 

-108 

0.139 

Flare 2 
56356 

32 

1.0 

1.3 

23.05 

-109 

0.157 

56355 

30 

2.7 

1.9 

14.60 

-78 

0.180 

56385 

29 

8.1 

2.7 

10.45 

-79 

0.157 

56384 

45 

5.8 

2.4 

8.35 

-82 

0.130 

56381 

65 

3.8 

2.2 

4.85 

-30 

0.132 

56433 

66 

10.2 

3.1 

2.64 

0 

0.086 

56463 

63 

24 

4.3 

32.00 

-73 

0.121 

Flare 3 
55747 

30 

1.0 

1.3 

14.60 

-45 

0.146 

55775 

29 

6.5 

2.5 

10.45 

-8 

0.114 

55812 

45 

4.2 

2.2 

8.35 

-29 

0.108 

55791 

64 

3.1 

2.0 

4.85 

-14 

0.108 

55806 

64 

8.9 

3.0 

2.64 

0 

0.084 

55820 

63 

24 

4.2 

32.00 

-24 

0.116 

Flare 4 
55638 

30 

1.0 

1.3 

14.60 

-26 

0.158 

55636 

29 

7.1 

2.6 

10.45 

-23 

0.145 

55639 

45 

5.4 

2.4 

8.35 

-25 

0.133 

55638 

65 

3.9 

2.2 

4.85 

-2 

0.110 

55661 

64 

9.0 

3.0 

2.64 

0 

0.081 

55662 

61 

25 

4.3 

32.00 

-73 

0.072 

Flare 5 
56183 

26 

0.8 

1.2 

14.60 

-57 

0.137 

56199 

27 

7.0 

2.6 

10.45 

-43 

0.128 

56213 

45 

4.7 

2.3 

8.35 

-15 

0.108 

56240 

65 

3.1 

2.0 

4.85 

5 

0.093 

56260 

64 

7.7 

2.8 

2.64 

0 

0.070 

56256 

61 

21 

4.1 

32.00 

-175 

0.400 

J0948+0022 

Flare 1 
55617 

44 

53 

7.6 

23.05 

-179 

0.320 

55614 

50 

69 

7.1 

14.60 

-153 

0.305 

55640 

61 

102 

8.1 

10.45 

-148 

0.249 

55644 

103 

57 

6.7 

8.35 

-92 

0.198 

55700 

103 

76 

6.6 

4.85 

-87 

0.098 

55705 

146 

51 

5.9 

2.64 

0 

0.064 

55792 

268 

34 

5.3 

32.00 

-168 

0.337 

Flare 2 
55705 

42 

49 

7.4 

23.05 

-159 

0.332 

55715 

50 

71 

7.2 

14.60 

-153 

0.309 

55721 

61 

103 

8.1 

10.45 

-128 

0.267 

55745 

103 

61 

6.9 

8.35 

-92 

0.185 

55782 

101 

73 

6.6 

4.85 

-65 

0.101 

55808 

146 

53 

6.0 

2.64 

0 

0.058 

55874 

260 

32 

5.2 
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Table B.2. continued. 


V 

(GHz) 

Time 

delay 

(d) 

Amplitude 

(Jy) 

Peak 

position 

(MJD) 

Time 

scale 

(d) 

T v ar 

(10 H K) 

D 

32.00 

-171 

0.331 

Flare 3 
56406 

42 

48 

7.4 

23.05 

-168 

0.277 

56409 

50 

60 

6.8 

14.60 

-133 

0.408 

56444 

62 

133 

8.7 

10.45 

-133 

0.400 

56444 

110 

81 

7.5 

8.35 

-125 

0.306 

56452 

108 

106 

7.3 

4.85 

-96 

0.138 

56481 

151 

68 

6.4 

32.00 

-173 

0.320 

Flare 4 
55446 

42 

46 

7.3 

23.05 

-161 

0.240 

55458 

50 

52 

6.5 

14.60 

-146 

0.119 

55473 

55 

48 

6.4 

10.45 

-139 

0.105 

55479 

93 

30 

5.5 

8.35 

-131 

0.094 

55487 

94 

43 

5.6 

4.85 

-96 

0.062 

55523 

137 

37 

5.4 

2.64 

0 

0.037 

55619 

230 

26 

4.9 

14.60 

-71 

0.092 

Flare 5 
54998 

52 

42 

6.1 

10.45 

-56 

0.071 

55013 

91 

21 

5.0 

8.35 

-41 

0.053 

55028 

87 

29 

5.0 

4.85 

7 

0.038 

55076 

128 

26 

4.9 

2.64 

0 

0.026 

55069 

206 

23 

4.7 

32.00 

-168 

0.180 

Flare 6 
55978 

37 

33 

6.5 

23.05 

-170 

0.163 

55976 

45 

43 

6.2 

14.60 

-175 

0.095 

55972 

52 

43 

6.2 

10.45 

-157 

0.098 

55989 

93 

28 

5.4 

8.35 

-147 

0.084 

55999 

92 

40 

5.5 

4.85 

-105 

0.059 

56041 

137 

35 

5.3 

2.64 

0 

0.032 

56146 

221 

24 

4.8 

32.00 

-104 

1505+0326 

Flare 1 

0.235 55545 

29 

48 

7.4 

14.60 

-37 

0.174 

55611 

50 

58 

7.6 

10.45 

-10 

0.177 

55639 

80 

44 

6.9 

8.35 

-13 

0.165 

55636 

78 

72 

6.9 

4.85 

-5 

0.148 

55643 

99 

111 

7.9 

2.64 

0 

0.112 

55649 

109 

233 

10.0 

32.00 

-120 

0.210 

Flare 2 
56146 

26 

52 

7.6 

14.60 

-33 

0.188 

56233 

50 

62 

7.8 

10.45 

-32 

0.175 

56234 

80 

44 

6.9 

8.35 

-30 

0.160 

56235 

78 

70 

6.8 

4.85 

-8 

0.156 

56257 

103 

108 

7.8 

2.64 

0 

0.128 

56266 

111 

258 

10.4 

32.00 

-107 

0.143 

Flare 3 
56088 

25 

39 

6.8 

14.60 

-55 

0.211 

56140 

50 

70 

8.1 

10.45 

-53 

0.197 

56141 

84 

45 

6.9 

8.35 

-49 

0.172 

56145 

78 

75 

7.0 

4.85 

-21 

0.138 

56173 

99 

103 

7.7 

2.64 

0 

0.101 

56194 

107 

219 

9.8 

14.60 

-27 

0.162 

Flare 4 
55809 

48 

58 

7.6 

10.45 

-29 

0.097 

55808 

77 

27 

5.8 
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Table B.2. continued. 


V 

(GHz) 

Time 

delay 

(d) 

Amplitude 

(Jy) 

Peak 

position 

(MJD) 

Time 

scale 

(d) 

T v ar 

(10 11 K) 

D 

8.35 

-16 

0.087 

55821 

73 

44 

5.9 

4.85 

-5 

0.092 

55831 

95 

75 

7.0 

2.64 

0 

0.078 

55836 

103 

181 

9.2 
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Fig.B.l. Left column: J0324+3410: top panel: the amplitude of the detected flares as a function of observing frequency with the 
cycle labels corresponding to the events seen in Fig.[2j middle panel: the flare delays; bottom panel: the 7\ al for the different activity 
cycles (flares) as a function of observing frequency. Middle column: J0948+0022. Right column: J1505+0326. 
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